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R-628 

AUTOMATIC LANDING SYSTEM OPTIMIZATION 
USING MODERN CONTROL THEORY AND 
INERTIAL MEASUREMENTS 


ABSTRACT 


Conventional automatic landing system designs for commercial transports 
have been primarily based on the utilization of Instrument Landing System (ILS) 
information to provide the position information required for acquisition and tracking 
of the landing glide path. Lateral position relative to a vertical plane passing 
through the runway centerline is derived from the localizer receiver output. 

Position relative to a vertical reference path is generated by processing the output 
of the gli deslope receiver. 

The spacial reference surfaces defined by the ILS signals are ideally plane, 
however, as a result of reflections from objects on the ground illuminated by the 
localizer and glideslope transmitter antennae, irregularities, referred to as "beam 
bending" or "beam noise", are superimposed on the planar surfaces. , Beam irregu- 
larities are the most significant source of noise in a conventional landing system. 

The amplitude of the noise is large enough to produce significant performance 
limitations on the landing system as a result of filtering and gain limitations which 
must be introduced due to saturation limits on vehicle effectors. 

Inertial navigator position and velocity is relatively free from the high 
frequency noise which plagues ILS information, however, erroneous.low frequency, 
gyro-drift induced, variations in position and velocity exist. The position error 
after a typical flight is measured in nautical miles, for example. This error would 
appear to preclude the application of inertial position for flight trajectory control 
during the terminal phases of flight. The possibility of combining the desirable low 
frequency characteristics of the ILS sensor with the broad bandwidth of the inertial 
information then suggests itself. By applying a systematic mathematical estimation 
technique to the signals from the ILS and inertial navigator it is possible to correct 
the inertially measured position and velocity. The result is position and velocity 
information which is relatively free from noise compared to the information provided 
by the ILS. 

In order to obtain maximum advantage from the noise free corrected inertial 
position and 'velocity data a complete review of the lateral and vertical position 
control systems was carried out. Primary attention was focused on improving the 
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magnitude characteristics of the open loop transfer functions of the vertical and 
lateral position regulators to increase bandwidth and minimize the effects of en- 
vironmental disturbances. The ability of the control system to follow a desired 
trajectory in space was increased by feedforward compensation. Maximum resist- 
ance to disturbances is achieved if the control effectors are always operated in an 
unsaturated condition. Optimization theory is applied to synthesize trajectories 
which the vehicle can follow (in the absence of disturbances) without saturating the 
effectors. This report develops the new control system structures and compares 
the performances of new and conventional automatic landing systems. 

Duncan MacKinnon 

January 1969 
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CHAPTER 1 


INTRODUCTION 


1. 1 Introduction 

Routine all-weather landing of transport aircraft is one of the outstanding 
unsolved problems of present-day technology. The importance of the problem and 
its implications both for passenger safety and for economic operation of airlines 
have been recognized for decades. Now that the jumbo-type transport and the 
supersonic transport are soon to be introduced into passenger service the problem 
of all-weather landing may be considered to be a critical factor limiting the full 
utilization of these aircraft. 

Over the past ten years, important strides have been made in the direction 
of the goal of all-weather landing. Of these, the most important has undoubtedly 
been the actual implementation of aircraft control systems which use Instrument 
Landing System (ILS) information and which have successfully performed thousands 
of automatic landings. An outstanding example is the system developed by the 
Blind Landing Experimental Unit (BLEU) of the Naval Aircraft Establishment, United 
Kingdom (see ref (1)).* It is highly probable that the basic concepts of control in 
what might be termed a classical automatic landing system such as the BLEU system 
will provide the conceptual framework for the systems that will eventually be adopted. 
It does not appear, however, that systems exist which will provide the invariable 
accuracy, reliability and independence from external disturbances, which will be 
mandatory for routine all-weather landing of transport aircraft. 

The point of view of the study described in this report is twofold. The study 
is concerned with the implications of inertial navigation technology for the all-weather 
landing problem. As a corollary, the study is concerned with the application of 
control theory, both classical and modern, for the effective utilization of the inertial 
equipment. The study is directed particularly towards the development of an all- 
weather landing system for the supersonic transport (SST). The SST presents a 
particularly challenging landing problem as a result of its high approach speed and 
its aerodynamic characteristics. It should be emphasized, however, that the re- 
sults of the study may be applied to any aircraft. 

^References (1) to (2 5) were referred to in reference 26 and are included in the table 
of references for completeness. 
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Over the past twenty years, an increasingly sophisticated technology of 
inertial navigation and guidance systems has grown up. This technology has been 
grounded on the continuous development of extremely precise reference components, 
principally floated gyros and accelerometers and a parallel development of support 
technology, such as gimbal structures, angular encoders, and thermal control sys- 
tems. The concurrent advances in computer technology have been indispensible to 
the rapid advances made in inertial systems. The developments in inertial navi- 
gation and guidance have been almost entirely motivated by military requirements. 
While the utility of inertial navigation systems for commercial aircraft has been 
recognized, the participation of commercial aircraft users in inertial navigation 
development has been a minor factor because of the costs involved. It is now con- 
sidered, however, that the accuracy, reliability and production costs of inertial 
navigation equipment originally developed for military pruposes are such that they 
are appropriate for commercial use. Plans for the large commercial transports 
soon to be introduced invariably include one or more inertial navigation systems 
per aircraft. It is planned, for example, that the SST will carry three inertial 
navigation systems to provide redundancy. 

Search of the available literature fails to show that any significant application 
of inertial navigation technology has been made to the automatic landing problem. 

A primary purpose of the present study, as noted, has been to demonstrate, by 
computer simulation, the implications of inertial navigation for the automatic landing 
problem. On the basis of work performed to date, it is believed that significant 
advantages will accrue from appropriate use of inertial navigation equipment. 

This contention is fully described and documented in the subsequent chapters. 

1. 2 Aircraft Landing 

The problem of landing an aircraft can be defined as a general problem in 
control theory. The control system, human or automated, must generate a set of 
commands which transfer the vehicle from some rather general initial state x^ to 
a terminal state x^. The initial state is usually on a linear trajectory inclined at 
some angle ifn relative to a vertical reference plane which contains the runway 
centerline. The terminal state, for the purposes of this discussion, is on a linear 
path described by the center of gravity of the aircraft in the vertical reference 
plane as it rolls down the runway centerline. The transfer between these states is 
usually divided into a number of clearly defined stages which are illustrated in 
Fig. 1. 2. 1. 
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Fig. 1.2.1 Nominal geometrical relationships at the initiation of localizer acquisition. 



1. Acquisition- During the acquisition phase the aircraft trajectory 

is changed to a path coincident with a reference defined by the 
terminal navigation system: In the case of an ILS system the 

path is specified by the intersection of two planes defined by 

the Localizer and Glideslope signals. 

2. Reference Path Tracking: During the path tracking phase the 
vehicle follows the reference defined by the terminal navigation 
system. The reference path is followed until the vehicle is 
70-100 feet above the runway elevation. During this phase a 
jet transport aircraft descends at approximately 10 feet/ second. 

3. Flareout: The flareout maneuver occurs during the final phase - 
of flight. The primary purpose of flare is to reduce the vertical 
rate of descent from approximately 10 feet/second to approximately 
2. 5 feet/second. Flareout is initiated at the termination of the 
Reference Path tracking phase. 

4. Decrab: If a component of wind perpendicular to the runway 

. centerline exists the aircraft is normally flown at a slight "crab' 1 
angle into the crosswind to reduce the lateral component of 
ground velocity to zero. The decrab maneuver, which occurs 
during flareout, aligns the longitudinal axis of the aircraft with 
the runway by rotating the vehicle in an uncoordinated fashion 
(roll <j> = 0). 

5. Rollout: Rollout Control is the final phase of vehicle guidance 

during a landing. The aircraft is on the ground, rolling along 
the runway. The rollout control system attempts to keep the 
center of gravity of the vehicle in the vertical reference plane. 

Thus the landing system must be capable of solving a relatively 
large assortment of control problems. The next section intro- 
duces criteria which may be applied to landing system performance. 

1. 3 Performance Criteria for Automatic Landing Systems 

Before the question of improving a system is approached, it is important 
to establish the exact nature of the indices which will be used to measure performance. 
These criteria should reflect the qualities which are normally desirable in flight 
control systems. In addition, special performance measures must be introduced 
which pertain to the problems peculiar to automatic landing. A set of measures 
which reflect these goals is: 
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Table 1.3-1 ' 


List of Specific Performance Requirements 


Phase 

Requirement 

Acquisition 

1. Acquire the ILS localizer and glide-slope centers 
. as quickly as possible with minimum overshoot. 

2. Perform this maneuver within the restrictions 
imposed on roll and roll rate. 

ILS Reference 
Line Tracking 

1. Minimize the error between the actual path of the 
aircraft and the ideal location of the ILS reference 
line. 

Flareout 

1. Minimize the error between a desired vertical 

velocity profile and the. actual vertical velocity profile. 

Decrab 

1. Minimize the lateral components of aircraft velocity 
and position at touchdown. 

2. Minimize the angular difference between a* vertical 
plane through the runway center line and the air- 
craft's longitudinal axis at touchdown. 

3. Attempt to achieve zero roll angle at touchdown. 

Rollout 

1. Minimize the distance between the path of the aircraft 
and the runway center line . 
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1. Sensitivity to environmental disturbances. 

2. Accuracy of flight relative to a desired reference trajectory. 

3. Control effector activity caused by noise. 

4. Physical limitations imposed by the aircraft structure. 

5. Human factors. 

The first two performance measures deal specifically with landing accuracy. 

A reduction in the effects of external disturbances ensures that landings may be 
repeated with small dispersion in spite of large changes in ambient atmospheric 
conditions. The ability to track a desired path accurately is obviously a particularly 
important item and will undoubtedly provide a key to winning the confidence of air- 
crews and the acceptance of the new automatic landing systems by the airlines. 
Accurate path tracking also provides the greatest margin of safety, since a com- 
paratively small deviation from the path may be interpreted as an incipient failure 
and an appropriate warning transmitted to the pilot. 

Control effector activity resulting from noise levels in the sensors which 
provide the information for control-loop closure must be restricted to a fairly low 
level to reduce wear on the effectors, decrease dragand limit undesired inputs into 
the pilots controls. 

The control system must also operate without exceeding the structural 
limitations of the vehicle. This implies control within a particular flight envelope 
and special care to ensure that the flexible bending modes of the vehicle are not 
excited. 

Human factors are particularly important in a passenger aircraft. Restric- 
tions on variables such as roll, roll rate and vertical acceleration must be incorpor- 
ated into the design. An automatic landing system should optimize items 1 through 
5 while working within these restrictions. 

The above criteria may be interpreted as a set of specific performance 
requirements applicable to each control stage of an automatic landing. Such a 
set of specific requirements is given in Table 1. 3-1. 

1. 4 Automatic Landing Systems 

As indicated above, automatic landing systems have been under development 
for some time. A conventional automatic landing system is shown in Fig. 1. 4. 1. 

The terminal navigation aid (ILS) provides position information x^ relative to the 
approach reference trajectory. This information is usually quite noisy as a result 
of fluctuations in the reference path due to random reflections from objects on the 
ground illuminated by the localizer and glideslope antennae. The output from the 
ILS is fed to a coupler which operates on the signal to provide inputs m Q tothe auto- 
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Fig. 1.4. I Conventional automatic landing system. 









pilot which are linear combinations of the position error, the integral of the position 
error and the derivative of the position error. The output from the coupler m Q is fed 
to a signal processor which limits the signal to account for human factors and flight 
envelope restrictions. The control signal mg is then applied to the aircraft autopilot 
which modifies the vehicle attitude to effect the trajectory changes required for ILS 
reference path tracking. The aircraft is placed under landing system control when 
it is initially on a linear path inclined at some heading angle iji a relative to the run- 
way centerline. The control system automatically acquires the reference path and 
guides the' aircraft along it until control is tranferred to a human operator or the 
decrab and flareout control systems at an altitude of 70-100 feet relative to the 
runway. 

While the simplicity of the conventional system is very attractive it suffers 
from a number of disadvantages. 

1. The open loop gains are restricted due to filtering which must 
be used to reduce the high noise level in the ILS signals. 

2. The control system tends to operate in saturation much of the 
time due to the high ILS noise level and the linear relationship 
between lateral position error and signal processor input. As 

a result the control system sensitivity to environmental disturb- 
ances is increased. 

The goal of this study is to improve the overall performance of a landing system by 
utilizing information provided by the inertial navigator and the contributions of 
modern control theory. This object is achieved by the application of four concepts 

a) Inertial filtering of the ILS reference signals 

b) Inertial stabilization of the aircraft control system 

c) Generalized trajectory control 

d) Nonlinear trajectory generation 

The concepts a-d form the basis for a new system which satisfies the criteria out- 
lined in section 1. 3. 

A block diagram of the MIT/IL automatic landing system is shown in figure 
1. 4. 2. The command input to the system is the instruction I which selects the 
operating mode of the landing system from the phases described in section 1. 2. 

The nonlinear trajectory generator develops a new spacial reference trajectory 
which interacts with the terminal approach reference trajectory to satisfy the 
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instruction I. The reference path r produced by the nonlinear trajectory generator 
accounts for 

1. Saturation constraints imposed by effectors and/or human factors. 

2. The bandwidth limitations of the closed loop position control 
system. 

As a result the vehicle is able to follow the trajectory precisely without saturating 
the effectors or the limits imposed by human factors. The utilization of trajectory 
generation has another very important implication since it permits the designer to 
select the system gains without considering the affect of gain values on gross re- 
sponse characteristics, permitting in general the utilization of higher gains. The 
output from the trajectory generator consists of a set of reference position, velocity 
and accelerations signals r. 

The position, velocity and lagged acceleration of the vehicle are generated by 
an integrated sensor package consisting of an inertial navigator, an estimator, and 
the terminal navigation aid (ILS). The output x c of the integrated sensor is a 
corrected set of inertial system measurements x^ which provide relatively noise 
free vehicle state data. The velocity, position, and acceleration components of x c 
are compared with corresponding components of r to generate a feedback control m Q . 
By exercising simultaneous control in this fashion over position, velocity and 
acceleration the affect of feedback path dynamics is eliminated and a considerable 
improvement in trajectory following precision achieved. A detailed discussion of 
nonlinear trajectory generation and its application can be found in chapter 4. 

In addition to the control component generated by comparing desired and 
actual values of position, velocity and acceleration a feedforward component q is 
produced to eliminate steady state errors. For example, a signal, in the vertical 
channel is generated to provide the vehicle pitch angle required to maintain a 
desired rate of descent. 

The total control signal is then applied to a signal processor which plays the 
same role as the command signal processor in the conventional automatic landing 
system. However as a result of careful trajectory generation and low sensor noise 
levels the limits in the signal processor are normally not exceeded. As a result 
the control system primarily operates in an unsaturated state reducing the sensi- 
tivity to environmental disturbances as explained in chapter 4. 

Subsequent chapters will discuss various aspects of the design of the 
MIT/IL Landing system in more detail.- 
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CHAPTER 2 


INERTIAL - TERMINAL NAVIGATION SYSTEM INTEGRATION 


2. 1 Introduction 

Inertially measured position and velocity are characterized by very low 
noise levels at high frequencies in spite of the wide bandwidth of the measuring 
instrument. Asa result inertial position and velocity would appear to be ideal 
feedback variables for position control system synthesis. Unfortunately the 
position and velocity information provided by the inertial system are subject to 
low frequency gyro excited errors which increase in amplitude with time. The 
magnitude of the error in position roughly increases one nautical mile for each 
hour of operation in a system of the type currently considered for commercial 
aircraft operation. Asa result the errors at the termination of a typical flight 
will preclude the use of position and possibly velocity as feedback control variables. 

A terminal navigation system such as ILS on the other hand, provides re- 
latively accurate position information which is corrupted by severe high frequency 
noise. ILS position and derived velocity are thus rather unsuited for control 
applications. 

The possibility of combining the desirable features of the ILS and inertial 
information then suggests itself. Such an integrated sensor unit would combine the 
low frequency characteristics of the ILS with the wide noise free dynamic range of 
the inertial system. The equations of integrated sensor providing corrected velocity 
and position information were developed in reference 26. This chapter extends the 
technique to encompass acceleration correction and examines the properties of the 
solution in more detail, 

2. 2 Initial Inertial System Correction 

The modern approaches to navigation system design utilize the total infor- 
mation provided by all the navigation sensors such as DME, VOR, VORTAC, 

LORAN, etc. and combines them using a systematic estimation technique such as 
Kalman Filtering or integral square error minimization. Asa result it is possible 
to update the inertial system during flight and it will probably be possible to reduce 
the position error of the inertial system below' one nautical mile before the aircraft 
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enters the terminal control area. Additional initial corrections may be made 
based on the data provided by the ILS system. Such initial corrections play an 
important role during the. acquisition phase of an automatic landing as shown in 
chapter 6 of reference 26. 

The lateral position relative to the vertical reference plane is the most 
important variable during the initial phases of an automatic landing. The lateral 
position may be estimated by examining the geometrical properties of a radius 
drawn from the ILS Localizer transmitter to the vehicle. Let r be the length of 
the radius and a the angular deflection of the radius from the vertical reference 

y 

plane. The lateral distance y from vehicle to vertical reference plane is then 
given by 


y = a r 

y y 


(2. 2-1) 


, Assume that , o> is measured by the localizer receiver while r y is generated by a 
distance measuring device (DME) located at the localizer antenna. The angle and 
distance measurements are subject to errors e and e resoectivelv. The actual 
measured distance y m is then given by 


(2. 2-2) 

The first term on the right hand side is the actual distance y. The last term is very 
small relative to the sec 
y is then approximately 


y = a r +»e +re +. e e. 
m y y ry y ay ay ry 


small relative to the second and third terms and may be neglected. The error e in 


e « s e + r e 

y y ry y ay 


(2. 2-3) 


If e and e are uncorrelated, the mean square value of the error E(e ) may be 
ay ry y 

written 


E(e r ) s* a 1 E(e ) + r 2 E(e ) 
y y ry y ay ’ 


(2. 2-4) 


Assuming typical values at the initiation of acquisition 


E(e ff y ) = 1 x 10' 4 rad 2 

a 2 - 36 x 10 4 rad 2 


E(e ry >.= 4 x 10 6 ft 2 
r 2 a 36 'x 10 8 ft 2 


12 



it is apparent that the second term in equation (2. 2-4) will dominate Efe^ ) an d 

E{e ) « r 2 E(e ) = 36 * 10 4 ft 2 (2.2-5) 

y y 

Since the mean square value of the error in inertial position E(e . ) on entering 
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the terminal area will range from 10 to 10 ft depending on the duration of flight, 
updating techniques and so forth, it is advantageous to perform an initial inertial 
position correction based on ILS - DME information by setting inertial lateral 
position y. = y m . 

Velocity data provided by the ILS - DME terminal navigation system, on 
the other hand, is generally much less accurate than the inertial velocity infor- 
mation. Thus an initial velocity correction based on terminal navigator data would 
not normally be made. 

A - similar approach may be used to correct the vertical position if it is 
measured by the inertial system. In this case it is generally advantageous to use 
the radio altimeter and/or barimetric altitude as the reference variable since the 
errors in these measurements are generally much smaller than the errors in the 
glideslope signal. 

While instantaneous inertial velocity and acceleration correction using ILS 
and DME information is not feasible it is possible to generate very precise inertial 
system corrections using mathematical estimation theory. This approach is dis- 
cussed in the next section. 


2. 3 Estimation of Inertial System Errors 

While instantaneous correction of inertial system position in 2. 2 results in a 
reduction in position error a more sophisticated approach must be utilized if a 
further correction in inertial position and corrections in inertial velocity and 
acceleration are desired. Mathematical estimation theory provides the techniques 
for accomplishing these goals. A least integral squared error technique is a 
relatively simple approach to this problem which requires only a modest quantity 
of computer memory. 

Let a set of inertially measured variables corresponding to motion along 
one axis be represented by a vector 


x 


il 


x. 

i 


X 


i2 


X i3 


(2. 3-1) 
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where 


x. j is inertially measured position 

x.£ is inertially measured velocity 

x^ is inertially measured acceleration 


A simplified set of relationships which models the inertial system is then 


x il ' x i2 
x i2 = x i3 


(2. 3-2) 


subject to the initial conditions x.^0), x i2 (0). Since the time interval during an 

automatic landing is short compared to the dynamics of the inertial system the 

errors in position and velocity may be. attributed to errors in the initial values 

x. 1 {0) and x.„(0) and an accelerometer bias error in x.„{0) {as a result of stable 
ll 14 10 

member misalignment and deflections of the local vertical). Thus 

x^(0) = (position at t = 0) + (error m position at t = 0) 

x.„(0) s= (velocity at t = 0) + (error in velocity at t =0) 

Id 

x i 3 < 0 ) = (acceleration at t =0) + (bias error) 

Knowledge of the errors in initial position and velocity and the bias error 
in acceleration is sufficient to correct all subsequent inertial data during the landing 
maneuver. 

Consider the system shown in Fig. 2. 3. !.• The inertial system described in 
equation (2. 3-2) is modelled by two integrators. The input to the inertial system 
is x. r. which is the sum of tnetrue acceleration x„ and the accelerometer bias b. 

id d 

The output inertial position x^j is fed to a model of the ILS receiver. The true 
position is fed to the actual ILS receiver. The outputs of the ILS receiver and its 
model are multiplied by range to the ILS transmitter r to convert to measured distance 
and compared at a summing junction. It is apparent that the resultant error e will 
be zero for all t > 0 providing 

1. The ILS model and ILS system outputs coincide at t = 0 . 

2. The ILS model accurately represents the ILS receiver. 
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3 The errors in inertially measured position, 
velocity and the accelerometer bias are zero 

4. The noise input n is zero 

It is assumed that the first two requirements are satisfied. Satisfaction of the 
third shall now be investigated. 

Let the vector describing the inertial system x. be perturbed at t =0 by 
6x^(0). Since the system is linear the effect 6e of the perturbation 6x.(0) on the 
error e may be written 

6e(t) = h(t ) ' 6x.(0) (2.3-3) 

where h is a vector of responses. The element h^(t) is the step response meas- 
ured at e to a step in x.^ applied at t = 0. The new error is then 

e(t> = e(t) + 6 e(t) (2.3-4) 

The estimation of the best perturbation in 6x.(0) may be accomplished by solving 
the following optimization problem 


Problem : Find a vector of perturbations 6x.(0) which minimizes the 

integral performance index 

t 


. J ~ J £e(<x) + 6e(a) J da 
0 

Expanding fe(a)+ 6e(a)J using equation (2. 2-3) gives* 

J = f |"e(a) 2 + 2h(a)' 6x.(0) e(a) 

J o 

+ Sx.(o)' h(a). h(a) # 6x.(0)Jda 


(2. 3-5) 


(2. 3-6) 


The first variation with respect to 6x.(0) in the performance index J is given by 


^The notation x' denotes the transpose of the vector x. 
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5 J = fix. (0) \ f 2h(a-) e(») 

0 L 

+ 2h(o) h(o)' 6x^(0) J dor 


(2. 3-7) 


The first necessary condition of the calculus of variations requires 6J to vanish. 
This condition is satisfied if 


6x.(0) = -TC h(or) h(o ) ' dol ^ h (a) €{a)da (2. 3-8) 

0 0 


The validity of equation (2. 3-8) becomes apparent if the ILS sensor noise is assumed 
to be zero: The error e(ar) then reduces to 


6(0-) = h(o) ' [x.{0) - x(0) J 

Substituting into (2. 3-8) yields 

6x^(0) = “jjy h(or) h(ar) ' dcrj - ^ 

% 

[x.(0) - x(0)]J 


{2 3-9) 


-1 „t 


h(a) h(o ) ' dar 


or 


x.(0) + 6x..(0) = x(0) 


(2. 3-10) 


(2. 3-11) 


the desired solution. The success of the solution hinges on the nonsingularity of 
the matrix 

t 


H = J [h(o) h(ar) '] do 


(2. 3-12) 


If t is small., t = 6t, integral 2.3-12 may be approximated by 

H = [h{0) h(0)'] 6t 

Since fit is a scalar the determinant |HI at t = fit is 

|H| - |h(0) h(0) '| fit 

Since the lh(0) h(0) 'j is zero the matrix H approaches singularity as 


(2.3-13) 
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(2.3-14) 



t -» 0 . This explains the transient behavior of the solution shown in Fig. 2. 5. 1. 

As time increases |H| will exist and a solution for 6x^0) generated. 

The noise n associated with reference signal x r introduces disturbances 
into the computation which affect the accuracy of the solution. An interesting 
property of the estimation algorithm is its ability to attenuate the effects of ILS 
noise n as time progresses. This effect is a result of the averaging process in- 
herent in the minimization of the integral (2. 3-5). The effect of ILS noise is 
illustrated in section 2. 7. 

The new initial conditions x^(0) + 6x^(0) may be used to calculate a corrected 
set of inertial information x c (t) by using the relation 

x c (t) = x.(t) + T(t) 6x^0) (2. 3-15) 

where the matrix T(t) is defined 

1 

0 
0 

A block diagram of this estimation procedure is shown in Fig. 2. 3. 1. This tech- 
nique may be used to correct vertical or lateral position velocity and acceleration 
information using the ILS localizer and glideslope signals. 

2. 4 Simplified Estimation of the Inertial System Position and Velocity Errors 

The general estimation procedure presented in the previous section may be 
simplified if the corrections are restricted to position and velocity. This is merely 
a specialization of the method developed previously and may be arrived at by rede- 
fining the following variables 





Neglecting the accelerometer bias the perturbation fie(t) becomes 


(2. 4-1) 


(2. 4-2) 
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(2.4-3) . 


fie(t) = h(t) ' fix^O) 


Substituting (2. 4-1), (2. 4-2) and (2. 4-3) into (2. 3-5) and minimizing the integral 
gives 


fix..(0) 



h(ff) h(o-) ' da 


ft 


hfo) e(a) da- 


(2. 4-4) 


The effect of an acceleration bias x^(0) creates an effective nonstationarity in the 
estimated initial value errors -fix.. The effect of acceleration bias is illustrated 
in the next section. The corrected inertial velocity and position may be computed 
once 6x^(0) is known 

x i (t) = x.(t) 4- T(t) 6 (0) (2.4-5) 

l 

where 


and 


x c (t) 



T(t) 


'1 

.0 


t 

1 . 


(2. 4-6) 


(2. 4-7) 


The response characteristics of this estimation technique are illustrated in the 
next section. 

2. 5 Response Characteristics of an Integrated Sensor Unit 

The behaviour of the estimation algorithm in section 2. 4 may be demonstrated 
’by individually introducing errors in position velocity and acceleration. The results 
of a sequence of runs of this form are shown in Figs. 2. 5. 1, 2. 5. 2 and 2. 5. 3. It is 
evident that the estimation process must continue for some time before useful infor- 
mation may be derived from the updated inertial position and velocity. To prevent 
the premature application of the information the corrections are multiplied by a 
matrix of functions U(t) before being applied to the inertial information. Thus 

x (t) = x.(t) + T(t) U(t) 6x.(0) (2.5-1) 

^ 1 1 
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Fig. 2.5.1 Corrected inertial position, y , and velocity, y , in response 
to an initial error in inertial position. (y=y = 0, y.= 0, y^= 500 ft) 
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Fig. 2.5.3 Corrected inertial position, y , and velocity, y , in- response 

C n c 

to an initial. acceleration' bias. (y = y = 0, y =0. 322 ft/sec^, y.= 0. 322t ft/sec 
y. = 0.161 t 2 ft) ‘ 1 



The matrix U(t) is a diagonal matrix. The elements u„(t) on the diagonal of U(t) 
are currently of the form 


u..(t) - 
11 


0 

t - u 


t <t x 


. V - 1. . h <t < t 2 

2i li 

1 t > t„ 


(2. 5-2) 


The times t^. and t 2 - are constants with respect to the initiation of estimation (t = 0 
in preceding sections of chapter 2). Equation (2. 5-2) shows that the correction 
6Xik(0) is gradually applied over the time interval (t^* *21^' "^is procedure 
reduces control system disturbances resulting from estimator transient's. 

2. 6 The Nature of Sensor Noise 

The primary source of sensor noise once estimator convergence is achieved 
in an ILS based integrated sensor unit is ILS beam center fluctuation. Changes in 
the glideslope and localizer beam center planes occur as a result of reflections 
from objects on the ground illuminated by the localizer and glideslope antennas. 
While the effect is deterministic in nature the three dimensional complexity of the 
pattern currently precludes the application of an analytic solution to eliminate the 
effects of the fluctuations. * As a result the perturbations in beam center position 
are usually characterized by an extraneous stochastic input to the ILS receivers. 

The receiver generally possesses a low pass filter which,, for the purposes of our 
investigation, is depicted by a simple time lag with time constant T . A typical local- 
izer receiver output resulting from beam noise is shown in Fig. 2. 6. 1. The power 
spectrum of this signal is illustrated in Fig. 2, 6. 2. 

In order to investigate the effects of beam noise on integrated sensor per- 
formance fluctuations in the beam center were simulated by exponentially corre- 
lating the output of a gaussian random number generator using a simple first order 
filter. The resultant signal in Fig. 2. 6. 3 was added to the input of the localizer 
receiver simulation to produce the output shown in Fig. 2. 6. 4. Figure 2. 6. 4 may 
be compared to the beam noise illustrated in Fig. 2. 6. 1 to verify the validity of 
the simulated noise. The power spectrum of the simulated receiver output is com- 
pared to the actual power spectrum in Fig. 2. 6. 2. 

Sensor noise level is conveniently quantized by generating the mean square 
value of the noise using equation (2. 6-1). 


*By storing the pattern, for example. 
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Fig. 2.6.2 Simulated and actual localizer power spectra. 



E (x) = limit I 

' ’ T-+ OP T 


( 2 . 6 - 1 ) 


I' 


x 2 do- 


If the random variable x is ergodic the limit (2. 6-1) exists and is equal to the 
2 

statistical mean of x . As a result of the relationships discussed in section 2. 2 it 
is apparent that the mean square value Efe^) of the error in y is related to the 
localizer output a^, the distance to the localizer transmitter r^. and the error e r y 
in Ty by the formula 


E(e ) 

y 


- «y E( 6 ry> 


+ r E(e ) 

y Q'y 


(2. 6-2) 


While E(e r ) and E(e^) may be considered constant during an automatic approach 
a and r are subject to variation. Thus the expectation E(e ) is time varying and 

y y 

the process e is nonstationary. Since O’ is usually very small after acquisition 

y y 

it shall be assumed that 


•E(e ) 

y : 


r E(e ) 

y «y 


(2. 6-3) 


A similar expression may be derived for the mean square value' of the error in 
vertical position E(e z ). 


E<e z > «,r/ E (e az ) (2.6-4) 

where r is the distance to the glideslope antenna 

is the angular deviation from the ideal glideslope reference plane 

The nonstationarity of e and with respect to distance indicated in equations 
(2. 6 2) and (2. 6 3) creates difficulties when two systems are to be compared. 

Thus it is more convenient to compare Efe^y ) and E(e^ z ) which may be generated 
from the relationships. 


> « r y - 2 E(e y ) (2.6-5) 

E(e ^ » - > < 2 6 ' 6 > 


2. 7 Integrated Sensor Noise Levels 

One of the most significant improvements achieved by sensor integration is 
a significant reduction in noise level. Noise originating in the ILS in introduced 
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into the inertial corrected" system output as a result of the correction algorithm in 
section 2. 4. 

- Since the noise is introduced in a linear fashion the noise contents in the 
corrected inertial position will vary as a function of the distances from the ILS 
localizer and glideslope antennas. In order to obtain meaningful statistical infor- 
mation stationarity must be re-established by dividing the corrected position by 
antenna distance before averaging. The response of the integrated sensor to the 
noise signals shown in Figs. 2. 6. 4 is illustrated in Fig. 2. 7. 1 The 
normalized responses are shown. The resultant averages may be compared with 
the values computed using equations (2. 6-4) and (2. 6-5). A similar procedure 
may be used to compare the ILS velocity estimates (obtained by differentiating the 
ILS position signals) with the corrected inertial velocity estimate. The resultant 
averages are given in table 6. 2-1. 

2. 8 Summary 

A systematic procedure has been presented for correcting "inertial system 
errors by combining ILS and inertial data utilizing mathematical estimation pro- 
cedure. A general algorithm is presented which simultaneously estimates errors 
in inertial position velocity and acceleration. This result is specialized to yield 
a simplified procedure for estimating position and velocity errors. The transient 
behavior of the solution is investigated by generating responses to initial errors 
in position, velocity and acceleration in the absence of ILS noise. The effects of 
ILS noise on the Integrated Sensor are introduced and the resultant position and 
velocity errors are compared to the corresponding errors in the ILS signal by 
generating mean square statistical averages. The results show that the noise level 
of the integrated sensor package outputs is quite small compared to the ILS signal 
noise. 
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Corrected inertial position divided by localizer range 



Corrected inertial velocity divided by localizer range 


Fig. 2.7. 1 Response characteristics of the integrated sensor 
unit to localizer beam noise, (t^ = t 12 = 15 . 0 sec, t 21 = t 22 = 25 . 0 sec) 
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CHAPTER 3 


LINEAR POSITION CONTROL SYSTEM 
SYNTHESIS AND PERFORMANCE LIMITATIONS 


3. 1 Introduction 

This chapter presents the pertinent elements of the synthesis and optimization 
of a new set of trajectory control systems which utilize the superior information 
provided by the integrated sensor package. The chapter also explores limitations 
imposed on the vertical and lateral position regulators by the dynamic characteristics 
of the vehicle and its effectors, the structure of the feedback control laws and the 
noise characteristics of the sensors. Error between the desired and actual values 
of controlled variables, sensitivity to disturbances emanating from the environment 
and control surface noise levels provide the basic criteria for this investigation. 

3. 2 Inertial Stabilization 

The integrated sensor unit, described in chapter 2, produces lateral and 
vertical position, velocity and lagged acceleration components which are relatively 
free from noise compared to the ILS localizer and glideslope signals. In order to 
effectively utilize this improved information a redesign of the lateral and vertical 
position control systems was undertaken. As shown in reference 2S the fundamental 
landing problem may be defined as a guidance problem in two orthogonal planes 
corresponding to the localizer reference plane, and a vertical plane containing the 
runway centerline. The intersection of the two planes defines the ideal ILS refer- 
ence path. This approach separates the control problem into two distinct areas, 
lateral guidance relative to the vertical plane and vertical guidance with respect 
to the localizer reference plane. As a result of the decoupling inherent in the 
simplified vehicle equations of motion an accurate system design may be achieved 
by separating the two problems using the techniques presented in chapter 5 of ref- 
erence 26. A convenient set of coordinates is an earth fixed reference frame 
originating at the point where the glideslope intersects the runway centerline. The 
x-axis lies along the centerline, y is horizontal while z is directed downward. The 
fundamental linearized relationships on which control system synthesis is based 
are then 
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(3.2-1) 


where 


y « g-cp 

z ~ Vp (a- 0) 


y is the lateral acceleration 
g is the gravitational constant 
cp is the roll angle 
z is the vertical velocity 
Vp is the path velocity 
a is the angle of attack 
0 is the pitch angle 


(3. 2-2) 


Lateral vehicle control is effected by modifying the roll angle cp using the roll angle 
autopilot. ' ; Vertical control results from modification of 6 and a using the 

pitch control system and/or direct lift control devices, such as spoilers. The 
pertinent linearized forward path transfer functions are 


and 



(3. 2-3) 



(3. 2-4) 


where 


^d ’ ^d are c * es * rec * r °H and pitch angles respectively 



is the roll control system transfer function 



is the pitch control system transfer function at 
constant path velocity 



is a transfer function relating the perturbation in 
angle of attack a to 0^ 
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The control variables, cp^ and 0^ are constructed by linearly combining the outputs 
of the integrated sensor unit, a feedforward compensation signal and a desired 
reference trajectory. 

FEEDBACK COMPENSATION 


^d 


r K y (K iy + s) x, K y 

* + K y s + r L rT I y c 


„ 2 -i 

K** s 

y 

r +~T 
as J 


FEEDFORWARD COMPENSATION 


K (K. + s) 

y iy 


+ K» s + K*» s 

y y 



(3. 2-5) 


0 


d 


FEEDBACK COMPENSATION 


K (K. + s) 

z 1Z 


K«» s 2 

+ K z s + ifn 
as 


FEEDFORWARD COMPENSATION 


K (K- + s) 


Z 1Z 


4- K* S 4- K** s 
z z 


z 

c 


where 



(3. 2-6) 


y , z 

J c' c 


y j z 

J r r 


V K * 

K iy' K i. 

K.. K* 
y z 


K j?' k*. 


are corrected values of y and z 
are desired values of y and z 

are adjustable position control loop parameters 
are adjustable integral compensator gains 
are adjustable velocity control loop parameters 

are adjustable acceleration control loop parameters 
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The differentiations shown in the feedforward component are avoided 
by explicit computation of the required derivatives as shown in section 4. 3. The 
structure of the reference components is based on the analysis in section 4. 4. 

(The small time constant T is neglected. ) 

3 - 

Thus it is apparent that the control signals cp^ and 0^ are formed by simul- 
taneously comparing the desired and actual values of position, velocity and accel- 
eration. Any discrepancy generates a command which corrects the error. The 
integral compensators (K^/s, K iz /s) weigh errors in position heavily as time 
progresses reflecting the emphasis on precise position'control. The feedforward 

n 

component s y /g provides the roll angle required to maintain the desired lateral 
acceleration based on relation 3. 2-1. The constant component of pitch required to 
maintain a steady rate of descent is generated by the term sz r /Vp. 

Since the structure of the control laws is fixed, analysis and synthesis leads 
to the definition of the feedback component parameters. Thus the designer is pre- 
sented with a parameter optimization problem. The criteria and constraints associ- 
ated with this problem are the topics of the next two sections. 

3. 3 Disturbances, Errors, Bandwidth and Open Loop Gain 

Two important criteria discussed in section 1. 3 are minimization of the 
effect of environmental disturbances on the system and the ability to follow a tra- 
jectory precisely. These phenomena are investigated by studying linear relation- 
ships using Laplace transform theory. A typical control system is shown in Fig. 

3. 3. 1. The object-of the control system is accomplished if the controlled variable 
C follows the reference input R. The system is subject to disturbances D which 
interfere with the control'process. The fixed elements' in the control system are 
and G 2 while Gq and H are subject to modification. 

The transfer function between C and R is 


C G 0 G 1 G 2 

R = 1 +GQGJG2H 


(3. 3-1) 


The dynamic response characteristics of the system depend upon the location of 
the zeros of the denominator. The location of the zeros may be modified by altering 
the transfer functions Gq and H. The zeros are normally located to yield an over- 
damped step response. The rapidity of the response is primarily determined by the 
magnitude characteristics of GqG^G 2 H. . 

The ability to follow a trajectory precisely in space is an important require- 
ment for an automatic landing system. Control system accuracy is measured by 
constructing the error between the desired and actual responses. 
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FIXED ELEMENTS OF THE CONTROL SYSTEM 


Fig. 3. 3. 1 Block diagram of a control system wr 
controlled output, C(s) and a distu: 






E = R - C 


(3. 3-2) 


The transfer function relating E to R is given by 

E 1 

R 1 + CgG^GgH 


(3. 3-3; 


The effect of external disturbances must also play an important role in 
control system design. Such effects are considered by examining the transfer 
function between the controlled output C and the disturbance D. 


C _ G 2 

D 1+G 0 G 1 G 2 H 


(3. 3-4) 


The open loop transfer function of the system is defined by the expression 


G ol “ G 0 G 1 G 2 H 


(3. 3-5) 


One of the most important control system parameters is the Open Loop 
Gain which may be defined as follows: 


Definition: The Open Loop Gain of the control system shown 
in Fig. 3. 3. 1 at s = Sq is defined by the limit 


G 


ol 


limit 
s -» s 


0 


G 0 G 1 G 2 H 


(3. 3-6) 


Asa result of the magnitude signs the value of the limit is independent of the 
direction in which Sq is approached if G^ is a rational fraction. 

In most physical systems the response to a sinusoidal signal of frequency u 
is of particular interest. In this case the general gain definition in equation (3. 3-6) 
reduces to the more common form 



limit 
s -» jcu 



(3. 3-7) 


The effect of variations in the magnitude and phase of the open loop transfer 
function may be determined by evaluating the magnitudes of {3. 3-3) and (3. 3-4). 



1 + G 


ol 


(3. 3-8) 


I C I 1 °2 

IdI=I i + g o1 


(3. 3-9) 


34 



Three cases are of interest 


1. !G ol t»l 


2- fGol' 5 * 1 


3 I G q1 I « 1 


E. 


The value of | 1 + G q ^| m | G^-,1 and j p- 


ol 


R 1 


ol 


C| 

d ! 


q — | . An increase in G q j in this range 


reduces the error and sensitivity to disturbances. 
The variation in | 1 + G q -J depends on the arg ‘G , . 

| ||| and j ^-| may increase, decrease or remain 
stationary as| | varies. 

| 1 + G q ^| ml, |^| and.| g-| are essentially inde- 
pendent of variations in|G o ^j 


The first case is the most important since it defines the spectral range where 
beneficial effects are obtained by feedback control. It is apparent that the frequency 
range where [ G q j| > 1 , the bandwidth, should be as wide as possible. The) G Q ^| 
should also be as large as possible within the bandwidth of the system. Thus the 
designer is presented with two related optimization goals The maximization of 
bandwidth and open loop gain is subject to a number of constraints which are con- 
sidered in the next section. 


3. 4 Limitations on Open Loop Gain and Bandwidth 

In general it is not possible to increase the open loop gain indefinitely. The 
factors which limit the open loop gain are 

1. The basic dynamic characteristics of the aircraft, effectors, and the 
signals available for control system stabilization. 

For a given set of fixed elements and a set of signals 
available for stabilization, there is a limit on the 
open-loop gain that can be achieved if the system is 
to be stable. 

2. The permissible control effector activity resulting from sensor 
noise. 

Residual noise in the position and stabilization signals 
further restricts the gain because of the limitations 
that must be placed on effector activity due to effector 
saturation and permissible levels of control energy. 
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3. The accuracy of the linearized model. 


Approximations are generally involved in the deri- 
vation of a linear mathematical model. Higher 
frequency dynamics which are usually ignored in 
the construction of a model become increasingly 
significant as the open loop gain is increased 
possibly invalidating the results of analysis. 

The stability of the control system shown in Fig. 3. 3. 1 may be explored 
using the Nyquist Stability Criterion. For the purposes of our investigation it is 
sufficient to assume that the open loop transfer function is equal to unity at one 
discrete frequency u>q and the magnitude approaches zero as u approaches infinity. 
In this case the Nyquist criterion states that the control system is stable (i.e. , all 
the zeros of 1 + G q1 lie in the left half of the complex plane) if the arg G ol is less 
than 180° and arg G qX is negative at u = u^. 

Consider the transfer functions G^. of the fixed elements and G y of the vari- 
able elements in the control system shown in Fig. 3. 3. 1. 


Since G v = G^ 
provided that 


G f = G x G 2 


G v - G 0 H 


at a unique frequency it is apparent that the system 


IgJ > |G f | -1 u < 


|G 1 - iG, 


i-l 


arg G y > - arg G f - 180 

a 




arg G y G f < 0 


> co = 


|G V I < |G X I (J > Wq 


(3. 4-1) 
(3. 4-2) 


is stable 


(3. 4-3) 
(3. 4-4) 


(3. 4-5) 


While equations 3. 4-3 to 3. 4-5 guarantee stability the transient response of the 
control system will be unsatisfactory unless tne angle [180 - arg GyG^Jat u)q , the 
phase margin, has a positive value. A phase margin of 30° to 60° is generally 
required to insure satisfactory transient response characteristics. 

The second limitation pertains to the level of control effector activity re- 
sulting from random noise originating in the vehicle sensors as a result of internal 
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processes and stochastic environmental effects. The noise level is usually meas- 
ured by evaluating the time average of the value of the noise squared. 

T 

E(x) = ^ ^ x 2 dt (3. 4-6) 

If the random processes are ergodic the integral {3. 4-6) will exist and will equal th 

statistical expectation of the squared value of x. Since E(x) is proportional to the 
2 

average value of x rather than the deviation in x itself it is usually more convenier 
to compare the square root of E(x) which is designated e(x) and referred to as the 
root mean square value of x in this report. 

Let 6 be an effector deflection. As a result of effector saturation and 
control energy constraints the mean square value of effector deflection must be 
limited. 

E<6> <E(6) max (3.4-7) 

The mean square value of effector deflection may be expressed as a function of 
the mean square values associated with the sensors and the parameters associated 
with Gy . Thus in order to satisfy inequality (3. 4-7), restrictions must be impose 
on the magnitude of the parameters in G v . The character of these constraints will 
become apparent when specific control systems are studied in subsequent sections. 

The accuracy of a mathematical analysis is measured by comparing the 
actual behavior of the physical system with the predicted behavior based on the 
equations which are used to model the system. Since every physical system is 
distributed, it is impossible to exactly model it with a finite set of differential 
equations. However, the solution of a finite set of equations may be sufficiently 
close to the actual solution for a restricted set of inputs to provide a useful pre- 
diction of actual vehicle performance. Increasing the bandwidth of a control syster 
results in the excitation of vehicle dynamics (such as bending modes) which lead to 
significant differences between predicted and actual behavior. As a result a more 
refined mathematical model must be used if a valid analysis is required. 

3. 5 Lateral Control System -Synthesis and Performance Limitations 

A general synthesis procedure which defines limitations on the bandwidth 
Uq and the open loop gain is now illustrated by application to the synthesis of 

a lateral position control system. A control system may be separated into fixed an 
variable sections. In the case of the lateral control system the transfer function 
relating lateral displacement y to desired roll is assumed to be fixed and is 
designated G fy . T his transfer function is shown in Fig. 3. 5. 1. 
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Fig. 3.5. 1 Magnitude and phase characteristics for the lateral 
position control system. 
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DE (db) 




(3. 5-1) 



The second area contains all the elements which the designer can modify to improve 
performance within the constraints discussed in section 3. 4. In this case the trans- 
fer function relating to y is subject to variation and is designated G 




s + 


K. 

iy 


+ K.s 

y 


+ ,_ K y 5 

T s + 1 
a 


{3. 5-2) 


v7hile the structure of equation (3. 5-2) is assumed to be fixed variation is made 
possible' by changing the values of the parameters K , K* , K.« and K. . 

y y y iy 

In the case of an automatic landing system the lateral position response to 
atmospheric disturbances is of prime importance. Response to disturbances and 
system stability is determined by the feedback component of the control law (3. 2-5). 
Therefore attention shall now be focussed on the definition of the adjustable feed- 
back parameters. A subsequent section will discuss the specification of a suitable 
reference trajectory. 

Consider the limitation imposed on bandwidth by stability. The lateral con- 
trol system is stable if the criteria imposed by equations (3. 4-3) to (3. 4-5) are 
satisfied. These criterias are satisfied if 


arg G vy > - arg G fy -180 


U s w 


oy 


(3. 5-3) 


The value of u,, is determined by the intersection of I G I and I G,. !” . Bandwidth 
Oy vy 1 1 fy‘ 

maximization is achieved by moving as far to the right as possible by manip- 


ulating alg G, 


vy 


It is apparent from the structure of (3. 5-2) that the bandwidth is 


maximized if the break frequency to^ y associated with the G vy is at zero frequency. 


This would imply that K y = 0, an unsatisfactory solution. 


Thus some compromise 


is required between maximizing t)Q y while maintaining a reasonable value for . 
A the following relationships exist between the 


Assuming that <s Wgy = 
gains K y , K> and K« : 


2y 


T 


!y 


1/2 


K. 


K. 


l k 5M 

K 

y ft 


2f, K 

iy y 
u iy 


(3. 5-4) 


(3. 5-5), 


(3. 5-6) 


(3. 5-7) 
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where £ ^ is the damping coefficient associated with the second order break fre~ 
quency u ly . A reasonable compromise is achieved by setting u ly = 0. 2 rad/sec. 
The .value of u Qy is then approximately 2. 0 rad/sec. The corresponding gains for 
t , = 1.-00 are 

b !y 


K 


y 


1. 13 degrees/ft 


K. 



11.3 degrees/ft 

2 

28. 3 degrees/ft/sec 


(3.5-8) 


The response characteristics of a lateral control system using these gains will be 

unsatisfactory as a result of the oscillatory nature of the system. The values in 

(3. 5-8) do, however, provide a well defined upper limit on the parameter magnatudes. 

Sensor noise introduces a further limitation of open loop gain and bandwidth 
by placing direct restrictions on the size of the adjustable parameters K y , and 

K y • 

While the general solution for the exact value of mean square control move- 
ment is extremely complicated, requiring precise knowledge of the statistical 
characteristics of all the disturbances acting on a system and the internally gen- 
erated design purposes using the following theorem. 

Theorem: Let a^ c** 3 l ^ ^ be any set of variables possessing finite auto- 

correlation and cross -correlation functions. Then the 
mean square value of their sum exists and is bounded 


E 


n n 

( > 1 i ) s 2 , 1 i 


E(a.) 


(3. 5-9) 


A proof of this theorem is given in Appendix B.. 

Referring to chapter 5 of reference 26 and equation (3. 2-9) the mean square aileron 
deflection may be , bounded 


where 


E(6 a ) < 2 K y 2 E<y) + K. 2 E(y) + kJ E(y) 


) 


+ E(cp) | K 2 + 2K| E<$) 


(3. 5-10) 


E(x) is the mean square value of the variable x 
K ■ are control system parameters 
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A typical set of worst case noise parameters for the new lateral control system is 


E(y) = 

400. 0 ft 2 

E(y) = 

25. 0 ft 2 / sec 2 

E(y) = 

4. 0 ft 2 /sec^ 

E(cp) = 

0. 01 deg 2 

E<$) = 

2 

0. 01 deg /sec' 


Suppose that the maximum permissible mean square aileron deflection is 

E(S ) = 2500. 0 deg 2 

el 


and the parameters associated with the gains are 


K = 1. 0 

cp 

K. = 2. 0 sec 
CD 


iy 

= 

y 


0. 200 rad/sec 

1 . 00 


The corresponding maximum set of gains may then be computed using eqs. (3. 5 -5), 
(3. 5-7) and (3. 5-3). It is apparent that the terms in equation (3. 5-10) contributed 
by roll and roll rate are comparatively small. Thus attention is focussed on the 
first three terms which are simplified using equation (3. 5-5) and (3. 5-6) to eliminate 


K* and K« 

y y 


therefore 


E(6 ) < 10800 K 
a' y 


K 


max 


E(6 ) 

a max 

10800 


(3. 5-11) 


K ^ 0. 481 degrees/ft 
K^.sa 4. 81 degrees /ft/sec 

K«sa 12. 02 degree/ft/sec 2 (3.5-12) 

Since the gains in (3. 5-12) are less than the gains in (3. 5-8) limitations on 
bandwidth and open loop gain are imposed by noise in this example. 


The integral compensator gain K.^ may now be selected. Since the integral 
compensator always introduces phase lag the break frequency u = K.^, must be 
selected so that the phase shift and magnitude effects introduced by the compensator 
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at the crossover frequency are small. This condition is satisfied if 


K. 

iy 


< ^ y 

10 


Thus the value of K. should be less than 0. 20. 

The final gains in table 5. 2-1 are based on the values in {3. 5-12) which 
were modified during simulation to achieve satisfactory transient response charac- 
teristics. The actual control system configuration is described in chapter 5. 


3. 6 Vertical Control System Synthesis and Optimization 


Conventional vertical position control systems utilize vehicle pitch to effect 

changes in the vertical plane. Thus the transfer function of the fixed elements G^ z 

in the control system relates vertical position z to desired pitch angle 0^ (eq 3. 2 -4). 

This transfer function is delineated in Fig. 3. 6. 1. The bandwidth limitations imposed 

by stability and noise are explored by plotting [G vz | 1 , - (arg G vz +180) and the 

control lay characteristics G imposed by the structure in equation (3. 2-6). This 

vz 

leads to the following conclusions 


1. The maximum crossover frequency Uq z lies between 1. 0 and 10. 0 rad / 
sec as a result of the phase restrictions on G vz - 

2. The slope of I G^ z I is approximately 40 db/decade' for 1. 0< u<10. 0. 
Since the maximum slope of G vz is 40 db/decade any increase in (j Qz 
in this frequency range, by acceleration feedback occurs at essentially 
constant open loop gain |G vz G^ 1 = 1.0. 

3. The break frequency w. should be as high as possible to maximize the 

1Z 

gain K z . 


Asa result of 2 the application of acceleration feedback to increase bandwidth will 
not lead to a significant decrease in sensitivity to environmental disturbances Thus 
acceleration feedback may be eliminated and the control law reduced to 


G = — = 
vz z 


K 


s + K. 

1Z 


+ K* s 
z 


(3. 6-1) 


Assuming that the integral compensator has a negligible effect at u 0z the maximum 

value of cj~ is approximately 2 rad/sec. The value of u. may be computed from 
UZ 

the relation 


w. 

1Z 




(3. 6-2) 


The integral compensator must produce negligible phase shift at Uq z - 
dition is satisfied if 


u 


K. < 

1Z 


Oz 

10 


This con- 
(3. 6-3) 
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ARGUMENT (DEGREES) 



LOG FREQUENCY 
RAD/SEC 


Fig. 3.6. 1 Magnitude and phase characteristics for the vertical 
position control system. 
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The upper bound on mean square elevator effector activity is estimated 
using the relationship 

E(6 e ) < 2 [K z 2 E(z) + K?' E(z)] K 2 

2 [K 0 2 E(fl) + K^ 2 E(fl) ] (3. 6-4) 

where K^,K^ are gains associated with the pitch autopilot described in reference 26 
and E(x 2 ) is the mean square value of the variable x. The second term in equation 
(3. 6-4) is comparatively small and may be neglected. 

Substituting 

O 

E(5 ) _ 900 degrees 

e max & 

E(z) = 100 ft 2 

E(z) = 25 ft 2 / sec 2 

K. = 0. 2 K (from equation 3.6-2) 

z z 

= 7. 10 

yields 

K z < 0. 351 degrees/ft 

This value of K is much larger than the level imposed by stability considerations 
z 

(K < 0. 287 degrees /ft). Thus the control system is stability limited and a value 

of K < 0. 287 must be used, 
z 

As a result of the limited performance improvement attainable using the 
conventional vertical control configuration the conventional approach was abandoned 
and a new technique which utilizes direct lift control adopted. The results of this 
investigation are presented in section 3. 7. 

3. 7 Improving Vertical Control System Performance by Direct Lift Control 

While vertical position control is not as critical as lateral position control 
an improvement in vertical trajectory accuracy,, particularly during the final landing 
phases, leads to important reductions in longitudinal touchdown dispersion. Inves- 
tigations in the past have dealt with the synthesis of a vertical ’path controller which 
uses pitch angle as the primary control variable. The sequence of phenomena which 
occur when a change in vertical path is effected by a pitch attitude modification may 
be summarized as follows 

1. An upward deflection command is fed to the elevator. 

2. The lift produced by the horizontal stabilizer decreases as the elevator 
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deflects up producing a positive pitch moment about the y axis. 

Since the SST horizontal stabilizer provides a significant portion of 
the total lift the vehicle begins to sink. 

3. The positive pitch moment leads to an increase in pitch angle which 
increases the angle of attack of the wing^ increasing the lift. 

4. In response to the increase in lift the vehicle begins to climb. 

It is apparent that vertical control by means of pitch attitude is subject to two im- 
portant disadvantages. 


1. The required elevator deflection produces an undesirable control re- 
versal at the initiation of the vertical maneuver. (This is a particularly 
serious problem during a landing abort maneuver. ) 

2. The relatively slow response of the pitch control system compared to 
the response of the vehicle effectors severely limits the bandwidth of 
the vertical position control loop. 

The disadvantages described above may be circumvented by effecting vertical con- 
trol using the direct lift spoilers. The direct lift spoilers are aerodynamic devices 
mounted on the top surface of the vehicles wings which may be deflected to disturb 
the air flow over the w'ings. This disturbance results. in a change in the total lift 
force. By operating with the spoiler deflection at a fraction of their maximum, 
deflection during the approach, positive or negative increments in lift may be generated 
by appropriate spoiler motion. The effect of spoiler operation is revealed by exam- 
ining the linearized vehicle equations for angle of attack a and pitch rate q given on 
page 234 of reference 26. The transfer function relating a to direct lift spoiler <f ^ 
deflection is derived in Appendix A. If the spoiler effector is modelled by a simple 
time lag with a 0. 05 second time constant T gi _j the transfer function between z and 
the spoiler command <^ S( ^ C may be written 


G 


z 

<5 , 
sdc 


T sd s+1 



q = 0 



(3. 7-1) 


This transfer function is shown in Fig. 3. 7. 1. The bandwidth limitations imposed 
by stability and sensor noise are investigated using the techniques in section 3. 5. 
The control structure is defined by equation (3. 2-6). Referring to Fig. 3. 7. 1 it is 
apparent that 

1. The maximum possible bandwidth is approximately 20 rad /sec 

compared to u)q z ps 4 rad/sec for the conventional system. 

2. The break frequency should be as high as possible to maximize 

the position gain K but u. cannot be much larger than 1 rad/sec if 
z iz 


45 



ARGUMENT (DEGREES) 



LOG FREQUENCY 
RAD/SEC 

Fig. 3.7. 1 Magnitude and phase relationships for a vertical control 
system using direct lift spoilers. 
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critical phase considerations are to be satisfied. 


3, . The application of lagged acceleration feedback only produces a 
minor improvement in bandwidth and open loop gain (u < u 0z ). 

Asa result the acceleration gain K«* is assumed to be zero. The control law in 

z 

equation (3,2 6) then reduces to 


sdc 


s + K. 

K — ^ — + K.s 

z s z 


(3. 7-2) 


The integral compensator gain K. must be selected so that phase shift introduced 

1Z 

by the integral compensator is negligible at This condition is satisfied if 


K. < 
lz 


Oz 

10 . 0 


(3. 7-3) 


The breah frequency is calculated from the relation 

K 


U lz = K 


(3. 7-4) 


The mean square spoiler motion as’ a result of externally and internally 
generated stochastic disturbances may be bounded using theoren 1 in section 3. 5, 


e( 6 . ) < 2 E < z > + K ; 2 

sac - z z 


Assuming that break frequency w - z is 1 rad /sec Kh is equal to and 


Substituting 


1 E<z) + E(z) 


E(6 gdc ) < 2 K z | E(z) + E(z) | 


E(6 , ) = 100. 0 

sdc max 


(3. 7-5) 


(3. 7-6) 


gives 


E(z) 

E(z) 


= 100 ft 

= 25 ft 2 /sec 2 


K = K- < 0. 63 degrees/ft 

Z Z “ 


This yields approximately 30 degrees of phase margin and acceptable transient 
response characteristics. Thus the final gain selection in once again based on 
effector noise considerations. 

While the direct lift technique improves the bandwidth of a vertical position 
control system its application presents a number of practical difficulties. 
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1 . 


The spoilers produce high drag. Thus it is desirable to operate 
as close to the reference setting as possible 

2 . The variation in z that can be achieved using spoilers is limited by 
magnitude restrictions on spoiler deflection. 

As a result of these problems a hybrid solution which combines the satisfactory 
low frequency response characteristics of the vertical controller based on pitch 
with the desirable high frequency properties of direct lift control is adopted. A 
system based on this configuration is presented in section 5.3 

3 . 8 Summary 

The performance limitations imposed on a linear control systems by fixed 
vehicle and effector characteristics, structural properties of the control law and 
mean square effector activity have~been explored. The results of this investigation 
have been applied to the synthesis of lateral and vertical position control systems 
for the SST. - While the properties of the linear controllers have been studied in 
depth nonlinear characteristics which have a profound impact on die nature of the 
vehicles path must now be considered. 
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CHAPTER 4 


OPTIMIZING TRAJECTORY PERFORMANCE 


4. 1 Introduction 

In chapter 3 it was demonstrated that the application of corrected inertial 
position, velocity, lagged acceleration data and direct lift control permits a signi- 
ficant improvement in position control system bandwidth. Asa result the sensi- 
tivity of the regulated position trajectory to ambient atmospheric disturbances was 
reduced and the ability of the controller to precisely follow a desired trajectory 
was increased. Bandwidth cannot be increased indefinitely as a result of consider- 
ations such as stability, noise and system nonlinearities. The most important 
vehicle nonlinearity during landing is saturation. Saturation occurs in all the 
aerodynamic effectors which produce translational forces and rotational moments 
on the vehicle. Additional saturation limits are imposed on variables such as roll 
angle, roll rate and' vertical acceleration by human factors such as passenger 
comfort. The implications of saturation nonlinearites in control system design 
are now considered. 

4. 2 Trajectory Characteristics and Saturation 

The most common form of saturation occurs in the effectors which produce 
the moments and forces required to effect changes in vehicle state. All the aero- 
dynamic effectors of the SST have magnitude and rate limitations. Effector 
saturation is particularly serious, resulting in an essentially open-loop condition 
as long as it persists. The effect of saturation is studied by examining the transfer 
functions of a system linearized about its current operating point, Consider the 
block diagram in Fig. 3. 3. 1 and let G 1 represent the effector transfer function. 
Saturation in effector output magnitude or rate modifies G^ to 

G 1 = 0 (4.2-1) 

Asa result, the transfer functions relating C to R and C to D become: 

C/R = 0 (4.2-2) 

C/D = G 2 (4.2-3) 
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Thus reference control is momentarily lost and the sensitivity of the system to- 
disturbances is increased by the factor 1 +Gq Gg H. It is therefore desirable to 
eliminate or prevent saturation. 

Saturation may be excited by three sources 

1. Reference commands to the control system. 

2. External disturbances acting on the vehicle. 

3. The combination of reference commands and external 
disturbances. 

Saturation arising from disturbances is controlled by reducing the magnitude of 
the open-loop transfer function. Thus the magnitude of the environmental dis- 
turbances imposes further restrictions on open-loop gain. Saturation resulting 
from reference inputs is restricted by processing the reference signal. This 
aspect is now considered in detail. 

4. 3 Nonlinear Trajectory Generation 


The problem of input signal processing may be formulated within the 
Theory of Optimal Control. For example, consider the linearized dynamical 
system 



x = f x + f m 
x m 

(4. 3-1) 

where 

X 

is an n-dimensional state vector 


m 

is an m-dimensional input 


f 

is an n x n matrix 


X 

f 

is an n x m matrix 


m 

The state vector is subject to a set of limits of the form* 



J 

VII 

7T 

(4.3-2) 


13 =S L. 

(4. 3-3) 


L and L. are n-dimensional vectors. The absence of a limit on a particular 

XX 

element of x is signified by setting the corresponding element of L equal to 
The problem may be defined as follows: 


*The inequality sign signifies that x. S L . i = 1, 2, . . . , n 

1 XI 
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Problem* Definition : Find a control m and the corresponding trajectory x, which 

minimizes the integral performance index 


T 

J = \ x' Mx dt 

t' 

0 


(4.3-4) 


subieet to the nnnhn'innnrQic constraint 

x = f x +f m (4.3-5) 

x m 

and the hard constraints 


ix| < L 1 
jx| < L 2 


(4.3-6) 


where 


L 1 < L x 
L 2 < L x 

The matrix M is positive semi- definite. 


(4.3-7) 


The object of the optimization is to transfer the vehicle to the vicinity of a 
terminal state while satisfying the limits imposed by Eq (4.3-5) The inequalities 
(4.3-6) are introduced to allow for the effects of disturbances on the system by 
providing some range between undisturbed operating values and saturation con- 
straints. 

The solution of this problem is complicated by the presence of the magnitude 
constraints. As a result, a two-point boundary problem must be formulated and 
solved. While the optimal control approach produces the best answer, the resulting 
computational complexity usually leads to a solution which requires a special- 
purpose, hybrid computer if real-time control is desired. These problems have 
discouraged the application of optimal- control techniques. 

To circumvent the computational difficulties the following approximate 
technique for generating solutions to the problem defined above is presented. The 
method is particularly suited to vehicle control problems. Consider the simplified 
vehicle control system shown in Fig. 4.3. 1. The desired terminal state may be 
set equal to 
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CONTROL PROCESSOR EFFECTOR 



s 


VEHICLE 


Fig. 4.3. 1 Illustrative model of a trajectory control system. The value of 
should follow the value of m^ closely. 










(4.3-8) 


x, = m 
1 c 

x 2 = 0 

The control input is the variable m^. The effector output x 2 is subject to 


4 2l 


< 




These limits imolv constraints on the vehicle state of the form 



Vi 


L 1 K 1 


Suppose that the natural frequency of the system in Fig. 4.3. 1 

j 1/2 
I *0*1 I 

w = 




is large compared to the bandwidth of the input m , and the damping ratio 



is greater than 1. Then the output satisfies 


provided that 



i* d | s b K i| 

s di 5 h K i 


(4.3-9) 

saturation 

(4.3-10) 

(4.3-11) 

(4.3-12) 

(4.3-13) 

(4.3-14) 

(4.3-15) 

(4.3-16) 

(4.3-17) 

(4.3-18) 
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The variable may be identified with the trajectory r. Thus the limits on the 
state variable x imply corresponding constraints on r if saturation within the 
control loop is to be avoided. A requirement is consequently established for a 
device which 


1. Controls the bandwidth of the input signal 

2, Modifies the input signal by constraining the maximum 
amplitudes of its first and second derivatives. 

A device having the desired properties is illustrated in Fig. 4.3. 2. The device 
will be referred to as a nonlinear trajectory generator. A linearized model of 
the NTG is obtained by removing the magnitude limits. The transfer function 
then becomes 


= K 1 
m c s 2 +K 2 K 1 s+K 1 


(4.3-19) 


The bandwidth of the transfer function is modified by varying the natural frequency 
and is normally equal to the crossover frequency of the open-loop transfer func- 
tion of the control system. The natural frequency depends upon the value of 

u t =K} /2 (4.3-20) 


The damping ration is normally larger than 1 


or 


? t ' 


tc K 1/2 
K 2 K 1 


2 ? 


(4.3-21) 


(4.3-22) 


The limits impose the following constraints on the outputs of the NTG 


l m dl= 

s o 

(4.3-23) 

N = 

S 1 

(4.3-24) 

N s 

S 2 

(4.3-25) 


Saturation of m^ automatically sets m^ equal to zero; similarly, saturation of m^ 
sets m ^ to zero. Thus the requirements imposed by saturation are satisfied. 
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The availability of first and second derivatives of the output plays an important 
role in the construction of control systems capable of precisely following a tra- 
jectory, as is shown later in the sections on vertical and lateral control system 
synthesis. 


4. 4 Generalized Trajectory Control 

The effectiveness of the approximate solution described in section 4. 3, 
depends upon the ability of the controlled system to precisely follow the generated 
trajectory. Trajectory precision has already been improved by 

1. Increasing the bandwidth of the control system. 

2. Controlling the bandwidth, damping ratio and saturation limits 
of the Nonlinear Trajectory Generator. 


The closed-loop transfer function between the reference input R and the controlled 
output of Fig. 3.2. 1 was stated earlier to be 


C_ G 0 G 1 G 2 
R 1 + G q G 1 G 2 H 


(4.4-1) 


It is of interest to consider the asymptotic behavior of this transfer function as 
the magnitude of the open-loop transfer function varies. Two cases are of interest 


(a) GqGjGjjH « 1 
R ■* G 0 G 1 G 2 


(4.4-2) 


(b) 


g q g 1 g 2 h » 1 



1 
H • 


(4.4-3) 


In the first case it is apparent that the advantages of feedback are lost if the 
open-loop magnitude is significantly less than one. In the second case the output 
C will not equal the input R unless H is identically one. This problem may be 
avoided by operating on the input with H so that Eq (4.4-1) becomes 


£ GgG^GgH 

R = 1 +G q G 1 G 2 H 


(4.4-4) 
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Consequently, precision trajectory control is achieved by producing a modified 
control signal which is a linear transformation of the desired trajectory. Typically, 
H has the form 


H — 1 + K* s + K**s 

X X 


2 


(4.4-5) 


where and Km are constants. The form of H in the present application implies 
that the first and second derivatives of the reference signal must be generated to 
provide the correct compensation. However, differentiation of the input can be 
avoided by utilizing the signals available from the nonlinear trajectory generator. 
The signals from the nonlinear trajectory generator are shown in Fig. 4. 3. 2 and 
this utilization is indicated in Fig. 4.4-1. The input R c to the control system has 
the required form 


R c = R d + K x R d +K x R d 

“ f 1 +K i s + K y s 2 )R c 


(4.4-6) 


- HR d 
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Fig. 4.4, 1 Block diagram of a cont 

trajectoi 










CHAPTER 5 


LATERAL AND VERTICAL CONTROL SYSTEM CONFIGURATIONS 


5. 1 Introduction 

Chapters 3 and 4 have discussed the theoretical considerations associated 
with inertially stabilized position control system analysis and synthesis. The task 
of the designer is to amalgamate these various concepts in a final design configuration. 
This chapter presents up to date modernized control system block diagrams and 
response characteristics. 


5, 2 Lateral Position Control System 

Lateral control is effective while the vehicle is acquiring and tracking the 
localizer center plane. During the decrab maneuver lateral position control is 
abandoned and lateral control responsibility is shifted to the Decrab Controller. * 
Lateral control is again resumed by the Rollout Control System* after touchdown. 
Rollout control has not been investigated to date. 

Lateral position control is achieved by turning the vehicle in a coordinated 
fashion (0 r* 0) using roll angle cp as the primary control variable based on the re- 
lationships 

• g tancp ( 5 2-1) 

' v 

P 


y « v p since 


(5„,2-2) 


Since the maximum amplitudes of cp and 4» are normally less than 45 degrees during 
the acquisition maneuver the response characteristics of the lateral control system 
and the linearized model shown in Fig. 5. 2. 1 are quite similar. The properties of 
the roll autopilot are described in appendix D. The parameter values associated 
with the control law are given in table 5. 2-1. Bode diagrams of the open and closed 
loop position controllers are shown in Figs. 5. 2. 2 and 5. 2. 3. One characteristic 
problem associated with the lateral control system is evident in Fig. 5. 2. 2. Since 


*See reference 26. 
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Table 5. 2-1 Lateral Control System Parameters 


GAINS 

K y 

Lateral position gain 

0. 080 deg/ft 

K y 

Lateral velocity gain 

0.800 deg/ft/sec 

K-. 

y 

Lateral acceleration gain 

2.000 deg/ft/sec 2 

K. 

iy 

Integral compensator gain 

0. 050 rad/sec 

CONSTANTS 

T a 

Lateral acceleration filter 
time constant 

0.100 seconds 

SATURATION LIMITS 

IMPOSED BY 
VEHICLE LIMITS 

IMPOSED BY 
NTG- OR CSP 

y 

CO 

CO 

♦ 

y 

±244 000 ft/sec 

±244. 000 


± 18. 60 ft/sec 2 

± 8. 000 ft/sec 2 

* 

± 30. 000 degrees 

± 30. 000 degrees 

* 

± 43. 790 deg/ sec 

± 10. 000 deg/sec 

y 

2 

± 21. 890 deg/sec 

2 

± 10. 000 deg/sec 

NTG PARAMETERS 

$ I 

aq 

Acquisition roll limit 

20. 000 deg 

aq 

Acquisition roll rate limit 

7. 500 deg/sec 

y 

y aq 

Acquisition roll acceleration 
limit 

2 

7 500 deg/sec 

INTEGRAL COMPENSATOR LIMITS 

L iy 

y integral compensator limit 

2. 0 degrees 
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Fig. 5. 2. 1 Linearized model of lateral control system. 
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the phase angle associated with the open loop transfer function approaches -270 
degrees as to approaches zero the lateral position control system is CONDITIONALLY 
STABLE. Control system stability depends upon the location of the crossover fre- 
quency . The value of w Q depends upon the gains K , K« , K-*, K. and the 

uy ny y y y y 

magnitude characteristics of the roll control system. If is sufficiently small 
the lateral control system will be unstable. While it is possible to adjust the control 
system parameters to obtain system stability saturation and trigonometric nonline- 
arities may effectively reduce the open loop gain to a level where a substained 
oscillation can occur. This problem is circumvented by limiting the output of the 
integrator in the proportional plus integral compensator so that the integrator will be 
in hard saturation whenever large values of roll are commanded. Integrator satu- 
ration effectively removes the integrator from the system reducing the phase lag to 
a safe maximum of 180 degrees at low frequencies. The step response of the line- 
arized system appears in Fig. 5.2.4. The effect of generalized trajectory control, 
discussed in section 4. 4, on closed loop bandwidth and step response is also shown. 
The wide bandwidth to reference inputs ensures that the lateral control system will 
be able to follow the desired trajectory precisely (in the absence of disturbances). 

The lateral position control system performs a sequence of distinct tasks 
during an automatic landing. The important features are identifed as 

1. Acquisition of the ILS Localizer center plane. 

2. Tracking the ILS Localizer center plane 

3. Momentary termination of lateral position control during the 
decrab maneuver. 

4. Resumption of position control during rollout. 

In order to accomplish these tasks the control system must transfer the aircraft 
from a linear path inclined with respect to the localizer reference plane to a path 
in the plane. The vehicle remains on this second path until rollout is completed. 

Lateral position control is subject to important limitations imposed by 
passenger comfort. These limitations apply to roll and roll rate. The trajectories 
generated by NTG^. must reflect the effects of these limitations. A convenient 
approach to lateral reference trajectory generation is based on the theory of time 
optimal control. A complete discussion of the synthesis of quasi -time -optimal tra- 
jectories is presented in chapter 6 of reference 26. The reference trajectories 
shown in Fig. 5. 2. 5 satisfy the relationships 

Kl * ^aq < 

I 'Pd * ®aq < 


c P 


max 


CD, 


max 


(5. 2-3) 


(5. 2-4) 
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1.0 


TIME 


30 sec 


Fig. 5.2. 4 Response of the linearized lateral position control system 
to step inputs in position error (solid lines) and desired position (dotted 
lines) with feedforward compensation. 
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Fig. 5.2.5 Acquisition control system responses for various initial 
heading angles. These responses demonstrate ability to acquire 
the localizer beam from a broad range of initial flight path headings. 


where cp and <p are the largest permissible roll and roll rate magnitudes 
respectively. and tp are constants. The absolute inequality in equations ' 

(5. 2-3) and (5. 2-4) provide a roll angle and angular rate margin for counteracting - 
the effects of environmental disturbances, tracking errors and inertial data correc- 
tions which are gradually applied during acquisition as shown in section 2. 5. The 
character of the trajectory and the wide bandwidth of the control system to reference 
inputs ensures that the aircraft follows the desired trajectory with a high degree of 
precision as shown in Fig. 5 2.6. Thus, if the environmental disturbances are 
moderate, the aircraft is able to fly the y reference path without exceeding the maxi- 
mum roll angle and roll rates. If the disturbances are very large, however, the 
commanded values cp c and <$> c may exceed the maximum permissible values <P max 

and m . In order to prevent excess roll and roll rate a nonlinear Command 
^max ^ 

Signal Processor CSP is utilized. The structure of CSP, is illustrated in Fig. 

cp 

4. 3. 2. The natural frequency of the CSP defined by equation (4. 3-20) and is adjusted 
so that linear response bandwidth of the CSP is much greater than the bandwidth of 
the roll control system. Asa result the output cp^ of CSP^ is essentially equal to 
the input c subject to the constraints 

Kl < Vax (5 - 2 ‘ 5> 


i*dl 


< 


q> 


max 


(5. 2-6) 


Thus, in absence of saturation, the CSP introduces neglibible effects on the linear 

response characteristics of the position control system It should be emphasized 

that CSP will normally operate in an unsaturated state thus ensuring minimum 
cp 

sensitivity to disturbances A comprehensive schematic diagram of the lateral posi- 
tion control system is shown in Fig. 5. 2. 7. 

5. 3 Vertical Position Control System 

Vertical control is maintained during all the phases of an automatic landing 
prior to touchdown. 

Vertical position control is achieved by modifying the vehicle pitch 6 and 
angle of attack or using the pitch angle control system and direct lift spoilers. The 
basic vertical plane relationship is 

z v sin{« - 9) (5.3-1) 

P 

The analysis of the vertical control system is based on the linearized model shown 
in Fig. 5. 3-1. " Since the perturbations in d and y are small during the landing pro- 
cedure the response characteristics of the actual control system are quite similar to 
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< 


Fig. 5.2.6 Typical acquisition control system response showing the reference 
trajectory y the actual path y and the roll angle <j> . = 45°). 
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Fig. 5 3. 1 Linearized model of a vertical position control loop. 








the linear characteristics. Magnitude and phase characteristics of the open and 
closed loop transfer functions are shown in Figs, 5. 3. 2 and 5. 3. 3, The step re- 
sponse of the linearized model appears in Fig. 5. 3. 4. The control system parameter 
values are given in table 5. 3-1. Pitch autopilot design is described in section 5. 11 
of reference 26. The important characteristics of the pitch regulator are illustrated 
in appendix D. The application of direct lift spoilers is illustrated in section 3. 7. 

As indicated in section 3. 2 the analysis of the vertical control system in based 
on the assumption that the perturbation in airspeed u is zero. Thus a critical part 
of the vertical control system is the automatic throttle. The automatic throttle 
utilizes airspeed information to generate an engine thrust command which maintains 
essentially constant airspeed thus satisfying the requirement {u « 0) for validity of 
the longitudinal transfer functions derived in appendix C of reference 26, The anal- 
ysis and design of an automatic throttle, of conventional configuration, is given in 
appendix E. 

The vertical trajectory may be divided into a three distinct segments 

1. Transition to the altitude h at which the localizer beam center 

cLCJ 

is intersected during the acquisition maneuver. 

2. Vertical flight path modification from a level trajectory h -h ac j 
to a descending path tangent to the glideslope centerplane. 

3. Reduction in vertical velocity from approximately 10 ft/sec to 
2, 5 ft /sec during flareout. 

This reference trajectory is shown in Fig. 5. 3. 5. In order to limit vertical accel- 
eration and position errors during transitions 1 and 3 particular care must be 
exercised in forming the trajectory. Let the vertical path be subject to the con- 
straints 


max 


. < z 
- max 


(5. 3-2) 


{5. 3-3) 


Assuming that precise control is exercised over vertical velocity and acceleration 
during acquisition and referring to Fig. 5. 3. 5 it is apparent that the following re- 
lationships hold for a smooth transition which satifies (5. 3-2) and (5. 3-3). 



< 


z 


max 


(5. 3-4) 
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' Table 5. 3-1 Vertical Control System Parameters 


GAINS 


' 

K z 

Vertical position gain 

0.063 deg/ft 

K. 

z 

Vertical velocity gain 

0. 126 deg/ft/sec 

K- ' 
z 

Vertical acceleration gain 

0.000 deg/ft/sec 2 
0. 050 rad/sec 

K. 

1Z 

Integral compensator gain 


50. 000 

K 

s 

Spoiler lead network gain 


CONSTANTS 



T 

Vertical acceleration filter 

0. 100 seconds 

ct •> 

time constant 



Spoiler lead network time 

5. 000 seconds 

s 

constant 


SATURATION LIMITS 

IMPOSED BY 
VEHICLE LIMITS 

IMPOSED BY NTG 

— 

z 

z 

— 

±14.600 ft/sec 

. z 

±25. 1 ft/sec 2 

±12. 500 ft/sec 2 

e 

— 

— 

• 

e 

± 5. 90 deg/sec 

— 

0 

± 4. 92 deg/sec 2 

— 

NTG PARAMETERS 




NTG z natural frequency 

0. 500 rad/sec 

^z 

NTG z damping ratio 

1. 000 


72 











ARGUMENT (DEGREES) 



MAGNITUDE (db) 





10 ° 10 1 10 2 


LOG FREQUENCY 
RAD/SEC 

Fig. 5.3. 3 Closed-loop frequency response characteristics to errors in position 
(solid lines) and reference inputs (dotted lines) with feedforward compensation. 




ro 

o_ 


1;0 


0 



TIME 


60 sec 


Fig. 5 . 3. 4 Response of the linearized vertical control system to step inputs 
in position error (solid lines) and reference position (dotted lines) with 
feedforward compensation. 
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T*'i g 5. 5 Nominal vertical trajectory. 



Ah b 
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1 

J 


(5. 3-5) 




2 


(x „ - x ) 
aq2 aq 


Ah 

a 

gs 


. 2 

x a 


X aql x aq2 


gs 


aq 


(5. 3-6) 


(5. 3-7) 


where 


« •* 

z , z 

max’ max 

are the maximum permissible vertical velocity 
and acceleration 

• • 

z aq 

is the desired vertical acceleration during 
acquisition 

V 

is the glideslope path inclination 

• 

X 

is the absolute vehicle velocity in the x 
direction. 

x , x x - ,Ah 

aq aq2 J aql 

are defined in Fig. 5. 3. 5 


The value of x & ^ - x a is approximately 700 ft. Thus the path transition 1 should 

be initiated approximately 3 seconds before the termination of acquisition at time 

T . The absolute inequality in equation (5. 3-4) provides an acceleration range for 
aq 

overcoming the effects of atmospheric disturbances, inherent tracking errors and 
vertical position corrections which are gradually applied as shown in section 2. 5. 
Reference trajectory vertical position velocity and acceleration serve as inputs in 
the complete control system shown in Fig. 5. 3, 7, As a result of the broad bandwidth 
of the vertical regulator to reference inputs precise tracking of the desired vertical 
path is assured. 

Initial deviations from the path described above are corrected using the 
nonlinear trajectory generator shown in Fig. 4. 3. 2. The errors in vertical position 
and velocity are introduced as initial conditions on z^ and z^ . The resultant outputs 
of the trajectory generator are added to the reference values z^ shown in Fig. 5. 3. 7. 
Thus the initial values of and z^ will correspond exactly to the estimated values 
of z and z. The response characteristics of the perturbation generator are shown 
in Fig. 5. 3. 6. 

The flareout trajectory may be synthesized in a similar fashion. As a re- 
sult of the importance of vertical velocity and acceleration regulation it is common 
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Fig. 5.3.7 Vertical position control system. 












to terminate altitude control during flareout. A detailed discussion of flareout may 
be found in chapter 7 of reference 26 

5. 4 Summary 

Block diagrams and performance data illustrating the current status of the 
lateral and vertical position control system designs have been presented. The 
performance of these position controllers will now be compared with the perform- 
ance of the conventional position control systems presented in appendix C. 



CHAPTER 6 


A PERFORMANCE COMPARISON OF CONVENTIONAL 
AND MODERNIZED POSITION CONTROL SYSTEMS 


6. 1 Introduction 

The significance of the innovations which have been introduced in the pre- 
ceding chapters can only be measured by comparing the performance of the mod- 
ernized control systems in chapter 5 with a set of conventional -control systems. 
Therefore an effort was directed towards the development of a conventional set of 
control systems for the SST which reflect the current state of the art in automatic 
landing system design. The conventional systems have been optimized using the 
same criteria applied to the inertially stabilized systems (with the exception of 
effector noise level criteria). The performance achieved by the conventional 
systems is representative of the capability of systems which provide path control 
utilizing position information derived from the ILS signal and velocity information 
based bn heading angle for lateral control and beam rate for vertical control. A 
discussion of conventional path controllers is found in appendix C. 

The preceding chapters have shown that significant improvements in per- 
formance can be expected in two -areas as a result of the application of integrated 
sensor information and modern control theory. These areas are 

1. Reduced sensor noise levels. 

2. Decreased sensitivity to environmental disturbances. 

The following sections will introduce criteria for measuring the performance in- 
crease. The application of these performance measures will demonstrate that 
significant improvements are indeed achieved by the application of modern tech- 
niques. 

6. 2 Sensitivity to Sensor Noise 

The integrated ILS Inertial Sensor described in section 2. 4 provides correctei 
position, velocity and lagged acceleration information. However as a result of the 
effect of the correction algorithm the velocity and position data is corrupted by ILS 
noise components superimposed on the corrected position and velocity. The noise 
components are nonstationary in the statistical sense as indicated in section 2. 6. 
Stationarity may be reestablished by dividing corrected position and velocity by the 
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appropriate distance to the terminal navigation system antennae. The sensor noise 

is quantized by comparing E («y ), E(y c / r y ) position and E(« y ), E (y c / r y ) 

for velocity while the vehicle is flown on a path y = y = 0 . The pertinent signals a , 

a , y /r and y fr are recorded in Fig. 2.7.1. The beam noise is shown in Fig. 
y -'c y^cy ° 

2. 6. 1. The mean square values are tabulated in table 6. 2. 1. These results show 
the dramatic reduction in noise amplitude achieved by inertial-ILS integration. 


Sensor noise in a conventional automatic landing .system is particularly 
serious as a result of its relatively large amplitude and the absence of any natural 
attenuating factor such as vehicle inertia which is very effective in the case of 
atmospheric disturbances. The importance of sensor noise may be inferred from 
the large perturbations which occur in the vehicle trajectories shown in Figs. 6. 2. 1 
and 6. 2. 2 as a result of beam noise. 


The effects of sensor noise may be quantized by evaluating the statistical 
averages of the noise generated errors. As a result of the unstationarity indicated 
above the signals' must be divided by the appropriate vehicle to ILS antenna 
distance. The values of E<— ) and E( Z ~ Z| 3 .) were generated during. simulation and 

r y r z 

are tabulated in Table 6.2-2. To facilitate the comparison the inertially stabilized 
variables are identified by the subscript "is" while the conventionally controlled 
quantities are subscripted "cs". 

The reduced sensor noise permits the application of higher feedback gains 
in the inertially stabilized control system loops compared to the gains in the con- 
ventionally stabilized system with similar response characteristics. Increased gains 
would tend to cancel the effects of sensor noise reduction, however, the important 
decrease in noise level still leads to a system which is comparatively free from 
beam noise generated trajectory errors as shown Figs. 6. 2. 1 and 6. 2. 2. ■ 


6. 3 Sensitivity to Random Atmospheric Disturbances 

The relationship between controlled output C and disturbance D of the con- 
trol system in Fig. 3. '3. 1 is given by 


C _ ... , G 2 (6.3-1) 

D 1 + 


Suppose that a comparison is desired between the responses of two different control 
systems characterized by variable element transfer functions of the form 


va 

" G 0a H a 

Vb 

= G 0b H b 


(6. 3. 2) 
(6. 3. 3) 
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Table 6. 2-1 

Mean Square Integrated Sensor Parameters 


e(o-y) deg 

0.224 

e (Ic)deg 

y 

0.100 

e(«y) deg 

0.224 

e (| £ ) deg 

y 

0.00511 



0.446 



0.022 
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The responses of the control systems to a disturbance D are given by 


G 2 

Ca = 1 + G 0a G l G 2 H a 

G 2 

Cb * TT V« 


The ratio of the responses is then 



1 + G 0b G l G 2 H b 
1 + G 0a G l G 2 H a 


I 1 


+ G 


J oll b 
ol a 


(6 3-4) 


(6. 3-5) 


(6. 3-6) 


Notice that the ratio of the responses is independent of the characteristics of the 
disturbance D. 

Suppose that C /.C, is plotted. Some general conclusions may be derived 
from the resultant curve. 

1. Jc^ / Ob 1 > 1 the response magnitude of system b is 
smaller than that of system a 

.2. j | = -1 the response magnitudes of the two 
systems are equal 

3. | G a/ G b[ < 1 the response magnitude of system a is 
less than that of system b 

As a result it is possible to compare the relative magnitude response characteristics 
of two systems. 

The statistical characteristics of the vehicle response to random disturbances 
play a key role in system performance valuation. If the spectral characteristics of 
the noise are assumed to be constant the ratios 


e < C a > 

e(D) ' e(D) 

and the relative performance criterion 


e^) 

e(< V 


(6. 3-7) 


(6 3-8) 
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Table 6. 2-2 
Mean Square Values 


Quantity 

Inertially Stabilized 
System 

Conventional 

System 

etoy) deg 

0.114 

0.114 

e(-Z) deg 

y 

0.038 

0.180 

e(a ) deg 
z 

0.105 

0.105 

e< j2L> deg 
z 

0.038 

0.120 


Derived quantities: 



b(%£)/ B(Oy) - 2.50 

e(U£)/e( Q z ) - 0.14 

! 2 » 

E(^a)/Eb z ) =1.30 
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•Conventionally stabilized lateral position 





Fig. 6.2. 1 Lateral control system response to localizer beam noise. 

^e(«y)l= 0.114 degrees j 
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I nertial ly stabilized vertical error 


e 


cs 


’IS 




Inertially stabilized system pitch angle 


Fig. 6.2.2 Vertical control system response to localizer beam noise. 

^e(« z ) = 0.121 degrees^ ' 
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may be evaluated. Quantities (6. 3-7) permit the general evaluation of response 
characteristics once e(D) is defined. The value of (6. 3-8) is very useful for com- 
paring different systems. The calculation of (6. 3-7) and (6. 3-8) is conveniently 
performed during simulation. 

The spectral performance of the lateral position control systems is compared 
in Fig. 6. 3. 1. A corresponding comparision for the vertical position regulators is 
shown in Fig. 6. 3. 2. The response of the vehicle to environmental noise (see chapter 
8 of reference 26) is illustrated in Figs. 6. 3. 3 to 6. 3. 6. The root mean square 
values of the controlled variables are summarized in tables 6. 3-1 and 6. 3-2. 

6. 4 The Effect of Windshear 

Windshear is defined as the derivative of the mean wind velocity with respect 
to altitude. Shear perpendicular to the runway centerline has a particularly adverse 
affect on the performance of the lateral position control system. The crosswind w 
is normally countered by the establishment of a crab'angle so that 

v sin s= - w (6. 4- 1) 

P y 

Since v); is small the following approximation is valid 


v p + - w y 


(6. 4-2) 


Windshear is perhaps the most serious disturbance encountered by the 
lateral position control system. In order to counteract the change in wind velocity 
with altitide, equation (6. 4-2) must be satisfied at all times. This implies that 


w 


4 » 


(6. 4-3) 


The heading rate is related to the roll angle cp by the relation 

• ct#uk 


gcp 

V 

p 


Thus the required roll angle is given by 

w 

* •* 

Since y and y are small during the localizer tracking phase 


cp =» K 


"s + K. "I 

*y 


y L s 


(6. 4-4) 


(6. 4-5) 


(6. 4-6) 


88 






25.0 ft 
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Conventionally stabilized position 



5.0 . 


ft/sec 



Conventionally stabilized velocity 



Fig. 6.3. 3 Lateral position control system responses to aerodynamic disturbance. 

^e{w n ) = 1.47 ft/ secj 
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Fig. 6.3. 4 Lateral position control system responses to aerodynamic disturbance. 

( e ^ w n ^ *1*47 ft/ sec^ 
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0 


I riertially stabilized vertical velocity 


41 secs 


Pig. 6. 3. 5 Vertical position control system responses in a turbulent .environment 

( ,e ^ w n^ = 1 • 46 ft/ sec j 
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Table 6. 3-1 


Root Mean Square Lateral Control System Parameter Values 
From Aerodynamic Noise Tests 


Quantity 

Inertially Stabilized 
System 

Conventional 

e(w^) ft/sec 

2.15 

2.15 

e(y) ft 

0.494 

1.11 

e(y) ft/ sec 

0,192 

0.297 

e(y)ft/sec^ 

0.132 - 

0.162 

e{6a)deg . 

0.607 

0.492 

** 


Derived quantities: 


e(y ) 
•^cs 


e (y is > 


2.25 


E(y. ) 
J 1S 

E(w n ) 


0.053 sec 


2 


E (y C s ) 

E(w n ) 


0.266 sec 


2 
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Table 6. 3-2 


Root Mean Square Vertical Control System Parameter Values 
From Aerodynamic Noise Tests 


Quantity 

Inertially Stabilized 
System 

e(w ) ft/sec 
n 

2. 22 

e(e ) ft 

Z 

3. 77 

e(z) ft/sec 

1.00 

e(6e) deg 

0. 656 

e ( 6 sd ) deg 

2.25 



2 . 22 


11. 81 


1. 70 


Derived quantities: 


e(e ) 
zcs 

e(e . ) 
zis 


E(e . ) 
zis 

E(w ) 
n 


= 2.87 sec* 


E(e ) 
zcs 

E(w ) 
n 


= 28.3 sec 
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multiplying both sides by s and transposing 


y 


S|£ 

J J iy 


(6. 4-7) 


It is apparent that the steady state value of y will be zero. The time required for 
the value of y to decay to 0.372 of its initial value is equal to K.y -1 ' The maximum 
perturbation in y due to windshear is bounded. An upper bound is given by 


H 


< l w yl 

max - gK~ 


(6. 4-8) 


if the response to a step change in wind acceleration is overdamped. Thus the 
performances of two lateral control systems to windshear are conveniently compared 
by evaluating {6. 4-8) and K. y _1 . The results of this comparision are contained in 
table 6. 4-1. The responses of conventionally and inertially stabilized systems are 
compared in Fig. 6. 4. 1. 

Windshear in the x direction produces errors in vertical position in 
similar fashion. If the perturbation in angle of attack a is assumed to be small 
and the airspeed is constant the flight path inclination may be approximated by 


v e 
P 

v + w 
P x 


(6. 4 9) 


The corresponding pitch angle 0 is given by 

y(v + w ) 

0 = P x (6.4-10) 

V P 

The path inclination angle yis determined prior to flareout by the constant elevation 
angle of the glideslope reference plane. Thus the value of 0 is constant in the absence 
of windshear. In the modernized control system the constant component of pis 
generated explicitly and applied as a control input as shown in equation 3.2-6. If the 
acceleration of the air mass in the x direction is not zero the pitch rate must obey 
the relation 

v w 

0 = ^ (6. 4-11) 

V P 

Suppose that w^ is constant. Then the pitch command must contain a component 
which changes linearly with time in order to satisfy equation (6. 4-11). In the 
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CD 




0 t 

I nertiaity stabilized lateral position 

Fig. 6. 4. 1 Vehicle responses to a step acceleration in air mass, of 0. 25 ft/sec 
which could result from windshear near the earth's surface. Note the different 'scales. 
The positive y perturbation is due to a minute initial lateral drift velocity. The 
steady state error is due to roll control system characteristics which are ignored 
in the linear analysis. 


absence of integral compensation the varying component of 0^ would be produced 
by a corresponding change in vertical error. With integral compensation the 
error is bounded and is given by the expression 


(z ' z d> = 


yw 


v K K. 

p Z 1Z 


(6. 4-12) 


The vertical windshear errors of the conventional and modernized control systems 
are conveniently compared by computing (z - z^/w^.) as shown in table 6. 4. 1. 

6. 5 Touchdown Dispersion 

A most important performance index, measures the error between the de- 
sired and actual touchdown points. Touchdown errors result from sensor errors 
(due to noise and inaccuracies) and aerodynamic disturbances. The variation in 
touchdown in conveniently investigated by considering the fluctuations in the lateral 
and longitudinal touchdown components independently. 

6. 6 Lateral Touchdown Dispersion 

The variation in lateral touchdown point .presents the simpler analytical 
problem. The errors may be divided into two components, a deviation due to beam 
and aerodynamic noise effects and a contribution due to windshear. If the various 
components are uncorrelated the mean square error may be expressed in the form 


E < e yt d ) = R(^. Oy) E(a y ).r ytd 2 + R(y, or)E(w n ) + R(y, w y ) E(£ y ) (6. 6-1) 

where e , , 
ytd 

r ytd 
R(u,v) 

w n is the stochasitc component of wind velocity. 

The first term on the right hand side of equation (6. 6-1} is the mean square value 
of the lateral position error resulting from localizer beam noise. The second 
term reflects the effect of turbulence in the atmospheric environment while the 
effects of windshear are introduced by the last term. 

A statistical windshear model was formulated assuming that the wind grad- 
ient is constant between h = 200 and h = 0 ft altitude. The mean square value of 
the deviation in y resulting from windshear may then be inferred from Fig. 6. 4-1. 


is the lateral position at touchdown 

is the distance from the localizer antenna at touchdown 
= E(u)/E(v) 
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conventionally stabilized 
system 

inertially stabilized system 

In order to present the preceding results in a clear fashion the value of E(w ) is 
assumed to be related to E(w ) by an expression of the form 

E(w ) = k* E(w ) (6. 6-2) 

y wy n 


R(y, wO 


9.0 4 

= 144.0 sec 


0. 0625 


1.0 ' 4 

= 16.0 sec 


0. 0625 


The value of k^ was set equal to 1.000. Equation (6. 6-1) may then be rewritten 


E V = 


r__x_, E(c? y )j 


+ R(y.j a ) + k^ ry R(y, w y ) j E(w n ) (6. 6-3) 


If it is assumed that the mean square value E(c- ) of the localizer noise is 0. 10 

2 ■ y 

degrees (reference 2) equation (6.6-3) may be plotted as a function of E(w^) 
using the coefficient values in tables 6.2-2 and 6.3-1. Curves illustrating the 
variation in Efe^^) as a function of ,w n are shown in Fig. 6.6. 1. 

6. 7 Longitudinal Touchdown Dispersion 

The analysis of longitudinal touchdown variation presents a more difficult 
problem. Touchdown occurs at the termination of the flareout maneuver which is 
generally initiated by sensing elevation above the runway using a radio altimeter. 
To simplify this investigation it shall be assumed that the error in longitudinal 
contact point is entirely due to the fluctuations in vertical path as a result glide- 
slope beam noise, aerodynamic disturbances and radio altimeter error. Assuming 
that the beam noise, atmospheric disturbances and altimeter error are uncorre- 
lated and the system response to the disturbances is linear the mean square value 
of the longitudinal contact error may be written in the form 


E(l 4d> -K(^.« z )E(Vf 2 «"Js 


[R(z,cr)E(w )+R . E(w )1 
L n zwx x j 


-2 


gs 


, -rr\f \ ■ - 2 2 

+ E{e hf )z td V P 


(6 . 7-1) 
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Table 6. 4-1 

Winds hear Comparisons 


Quantity 

Modernized 
. System 

Conventional 
System • 

(gKy) ^sec^ 

22.2 

87.0 

K iy _1 sec 

20.0 

200.0 

y 

V K K. sec 2 
p z 1Z 

0.873 

0.057 


Derived quantities: 

<S K y>i 



3.92 



100 




















where 


'xtd 


a 


gs 


w. 


'hf 


Td 


is the difference between the desired and actual touchdown 
points along the runway. 

is the distance from the glideslope antenna at the initiation 
of flareout. 

is the glidepath angle . 

is the stochastic component of wind velocity. 

is the radio altimeter error. 

is the nominal vertical velocity at touchdown. 


The first term on the right hand side of (6. 7-1) is the mean square error resulting 
from glideslope beam noise; the second term is the contribution due to atmospheric 
disturbances while the last term represents the contribution as a result of altimeter 
error. 

The properties of equation (6. 7-1) may be illustrated by setting 


a =0. 0437 radians 

gs 

r f = 1500 ft 

Eie 2 ^) = 4. 0 ft 2 

z w 2 = 6. 25 ft 2 /sec 2 


The above values may then be used to generate values of 6. 7-1 as a function of 
w n which are then summarized in Fig. 6. 7. 1. 


6. 8 Root Mean Square Touchdown Area 

The results presented in the preceding sections maybe applied to illustrate 
the touchdown performance of an automatic landing system in a graphic fashion. 
The root mean square touchdown area is defined by a width equal to 2e(e y1 . d ) and a 
length equal to 2e ( e 3C j ;c j)- These areas may be drawn using the plotted data in Figs. 
6, 6. 1 and 6. 7. 1 as shown in Fig. 6. 8. 1. 


Assuming that the deviations in touchdown point are normally distributed 
the following probability relationships may be written 
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: — i 100 ft 

conventionally stabilized 

jL 

■ 100 ft 

inertially stabilized j- 

e(W ) = 0.0 ft/sec 
n 


conventionally stabilized 

± 

1 100 ft 

inertially stabilized T 

e(W ) = 4.47 ft/sec 


Fig. 6.8.1 Root mean square touchdown areas for various values 
of mean square turbulence. 





p 


= 0 . 68 


( e ytd S e(e ytd )) 

P ( e ytd S 3 00 e <Vd ) ) = 96 

P ( e ytd ^ 2 ' 58e(e ytd ) ) =°- 99 (6.8-1) 

Thus, we are assured that an aircraft performing landings in a stationary (in the 
statistical sense) environment will land within the root mean square area 68% of 
the time. The maximum permissible deviations in y and x at touchdown are nom- 
inally defined as 50 ft and 1000 ft, respectively. In the example presented in sections 
6. 6 and 6. 7 lateral error is the critical parameter. The modernized control system 
performs 99% of the time with root mean square turbulence velocities in excess of 
10. 0 ft/sec. The conventional system.on the other hand,fails to achieve a satis- 
factory landing 32% of the time at a modest value of e(w^) =4.0. 

6. 9 Summary and Conclusions 

Thepreceding sections of chapter 6 provide the basis for some general 
conclusions concerning the development status of the NASA /MIT automatic landing 
system. The lateral control system will be discussed first followed by the vertical 
position regulator. 

If the. results in sections 6. 2, 6. 4 and 6. 6 are studied.it is apparent that 
the most important factors which lead to large lateral touchdown dispersions in a 
conventional lateral control system are 

1. High levels of ILS induced noise in the vehicle position 
and velocity data. 

2. Sensitivity to windshear. 

By integrating inertial navigator and ILS information a very significant 
reduction in measured position and velocity noise components was achieved. As a 
result the root mean square perturbations in lateral trajectory due to beam noise 
are reduced by a factor of 4.73 as shown in table 6. 2. 1. 

The lag-free nature and low noise characteristics of the corrected inertial 
lateral position and velocity data added to the availability of lagged lateral accel- 
eration information provided the basis for a new lateral position control system 
design which features important increases in the position feedback gain K^, the 
velocity gain K^., the integral compensator break frequency K.^ and a new loop 
closure on lagged inertial acceleration through a gain K" . By carefully tailoring 
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the desired trajectory of the vehicle to prevent effector and/or flight envelope 
saturation, maximum advantage can be derived from the increased open loop gains 
and response bandwidth to minimize the effects of atmospheric disturbances 
Satisfactory lateral reference traj'ectories are currently synthesized using the 
nonlinear trajectory generation techniques presented in chapter 6 of reference 26. 

The combination of increased gains and nonlinear trajectory generation resulted 
in an approximately 2 times reduction in the root mean square amplitude of the 
lateral path perturbations due to atmospheric turbulence. The most" important 
reduction occurred in the lateral perturbation at touchdown from wmdshear. The 
inertially stabilized system reduced the effect of windshear by a factor of 12. 

In section 6. 7 it was shown that the dispersion in the location -of the touch- 
down point along the runway centerline is primarily caused by 

1. Errors in the altitude at which the flareout maneuver is initiated. 

2. Deviations in the vertical trajectory as a result of atmospheric 
disturbances. 

A reduction in the first source of error is limited by the accuracy of the 
instrument (usually a radar altimeter) which provides the automatic landing system 
with runway surface referenced elevation data during the terminal phases of the 
landing. This error may be reduced by improved instrument accuracy and/or 
inertial system-altimeter integration using the techniques presented in chapter 2. 

The work to datehas not covered this facet of the landing problem. The effects of 
the second dispersion generating factor, atmospheric turbulence, may be reduced 
by increasing the open loop gain and bandwidth of the vertical position control 
system. The gain and bandwidth of the conventional system are limited by the 
noise levels in the measured vertical position and velocity data and the inherent 
physical characteristics of the vehicle. The first limitation is circumvented by 
combining inertial and ILS information. Sensor integration reduces the root mean 
square noise in position and velocity by 2 and 45 times respectively (table 2. 6-1). 

By applying auxiliary direct lift control the limitations imposed by vehicle char- 
acteristics are eased and the vertical control system open loop gain and bandwidth 
may be improved. The root mean square vertical deviations that result from stochastic 
components in the atmosphere are reduced by a factor of 3 by the improved open 
loop gam and bandwidth. 

The combined improvements in the vertical and lateral position control 
systems result in a landing system which is able to perform successfully m a 
much broader range of ambient turbulence and windshear as shown in section 6. 8. 
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CHAPTER 7 

AREAS FOR FURTHER INVESTIGATION 


7. 1 Areas for Further Investigation 

In the course of the program summarized in this report an attempt has been 
made to elucidate and solve some of the basic problems associated with automatic 
landing systems. The results presented in the preceding chapter indicate that 
some major improvements are feasible if inertial and terminal ITS navigation 
system data is properly combined and the control system design is revised to re- 
flect the improved character of this information In addition to these concrete 
results many significant problems have been revealed which could not be fully 
investigated as a result of time and funding limitations. Among the tasks particu- 
larly appropriate for future investigation are: 

1. Investigate the affect of inertial system quality on landing system 
.performance. Determine the lower limits on inertial system accu- 
racy for satisfactory landing system performance. 

2. Study the practical problems introduced by disparities between the 
locations of the vehicle center of gravity and the inertial navigator. 

3. Apply inertial -ILS-DME* filtering techniques to generate improved 
information along the glidepath. 

4. Utilize inertial system data to detect terminal navigation system 
failure and/or interference caused by overflight. 

5. Investigate the feasibility of extrapolating the reference glidepath 
using inertial .navigator information during periods of terminal 
navigation system failure and/or interference. 

6. Study the application of longitudinal direct force generators 
(dive brakes ) to improve the performance of the automatic 
throttle. 

7. Redesign the automatic throttle to incorporate inertial acceleration 
data. 

8. Improve the lateral control system by the application of side force 
generating effectors, 

#DME is not available at some facilities . 
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9. Optimize the landing abort procedure using variational techniques 
and design a landing abort control system. 

10. Refine the flareout control system design presented in reference 
26 to include the latest improvements in vertical control system 
structure and inertial- radio altimeter data processing. 

11. Synthesize an optimal "decrab" maneuver trajectory using mathe- 
matical optimization techniques and redesign the decrab control 
system. 

12. Model the vehicle in the rollout configuration and design a rollout 
control system. 

13. Investigate the application of inertial information to monitor and 
control aircraft takeoff and takeoff abort. 

14. Perform' flight tests to validate the results obtained using the 

, simulation. 

The availability of an on board inertial navigation system offers unique possi- 
bilities for flight testing new automatic landing system concepts. Landing system 
tests have historically been performed at low altitudes utilizing installed ILS facil- 
ities. As a result a certain degree of danger has existed in the initial test phases 
as a result of proximity to the ground. Since an inertial system provides an absolute 
geographic reference a hypothetical runway may be defined at an elevation above 
ground. The geometry of the runway, an instrument landing system, DME etc. may 
be defined. The location of the test vehicle relative to the imaginary airport is then 
defined by comparing the aircraft's inertial position to the geometrical coordinates of 
the runway complex. Computations may also be performed to generate the signals 
received from the simulated ILS and DME (including ILS and DME noise). Thus 
the aircraft could perform all the phases of an automatic landing at a safe altitude. 
The effects of turbulence are easily evaluated by locating the runway complex in an 
unsteady atmospheric environment. Simulation of the effects of wind or windshear 
is easily accomplished by translating the geographic coordinates of the hypothetical 
airport. .The absence of ground effect and 1 touchdown dynamics would contribute the 
primary source of error. Once the landing system is operating satisfactorily at 
high altitudes low altitude tests could be performed using actual ILS and DME 
information. 
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APPENDIX A 


SPOILER TRANSFER FUNCTION 


Derivation of the transfer function relating angle of attack to direct lift 
spoiler deflection follows 

The equations relating a and q to the elevon deflection 6 and the direct 
lift spoiler deflection 6 sd are given in reference 26. 

C.a = C. + C. q + C. .9 + C. a + C . 6 

a au fixp ae aa ax e 


+ C. * 6 . 

a6 sd sd 


(A. 1-1) 


where 


C.q 

q 


C. u 
qu 


+ C. qt C. a + C. , 

qq^ 1 qa qs. 


C q6 sd ^ 


(A. 1-2) 


o* 

a is a constant 

C • is a constant relating a to the variable x 


u is the perturbation in airspeed 

a, q are perturbations in the angle of attack and pitch rate respectively 
5 , 6 gd are the deflections of the elevons and direct lift spoilers respectively 

Assuming that u is zero as a result of the control effected' by the automatic throttle 
the equations reduce to 


C.a = C, q 
a aq 


C. 6 
aO 


C. a 
a a 


+ C 


a 6 


+ C. - 6 , 

a6 sd sd 


(A. 1-3) 


C.q = C. q 

q qq 


+ C. a 
qa 


-r C • c 5 + C. 


qfi 


«*sd sd 


(A. 1-4) 
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It is desirable to eliminate the pitching moment due to spoiler deflection. This 
goal is achieved if the elevon is deflected so that 


6 = - C._ " 1 [c. a + C., 6 j] 

e q6 e L qa q6 sd sdj 


(A. 1-5) 


As a result it is possible to set q and 9 equal to zero in equations (A. 1-1) and (A. 1-2). 
Substituting for 6 g in equation (A. 1-1) gives 

C.cc = [C. - C. t C. -1 C. 1 o 
a L aa a6 q6„ qaJ 


< 


+ C.. - C.„ C 


-1 


cc8 


a6„ q6 q8 c 


sd 


= C. oc + C. , 8 j 

aa ao gd sd 


(A. 1-6) 


Rewriting the equation in Laplace transform notation and rearranging terms yields 


[c.s - C. a = C. . 6 , 

L a aaj a6 gd sd 


(A. 1-7) 


or 


/ _a_\ _ ^aS S( j / C a a 

\ 6 sd/ C. /C. s - 1 
q=0 a 7 aa 


(A. 1-8) 


Using the stability derivatives 


c. 

a 

— 

- 4. 630 

C. 

aa 

= 

3. 235 

C a6 p 

e 

= 

0. 288 

C.* 
ao , 
sd 

= 

- 0. 731 

c ^« 

= 

- 0. 079 

c. 

qa 

= 

- 0. 043 

c.. 

^ 6 sd 

= 

0. 000 


and equation (A. 1-5) gives 



3. 080 


C. 

act 


C.. 

ao 


0. 731 


sd 


The transfer function (A. 1-7) is then 



0. 237 

1. 510s + 1. 000 


(A. 1-9) 
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APPENDIX B 


B. 1 Proof of Theorem 1 

Let the time average be defined 


i. 

A, 2. limit 1 P 2 
^ (X } = T->» T \ X dt 


(B.l)' 


and let the time averages of the squares and cross products of the elements a. 
i = 1, n exist. Consider the first two elements of the sequen.ee, 

■ ^(( a i + a 2^) = f?^ 2 ) + £(a 2 2 ) + 2^*{a^a 2 ) 

fi‘((a 1 - a 3 ) 2 W( ai ,) 2 +(?(a 2 2 ) T 

Since the left hand sides of B.2 and B.3 are positive semi- definite 
-<f( ai 2 ) +<?(a 2 2 ) >2(?{a 1 a 2 ) 

Therefore- 

<?(( ai +a 2 > 2 ) < 2^{ ai 2 ) +(?(a 2 2 )) 

Now suppose that the theorem holds for the first k terms (k<n) so that 




Then for k + 1 terms 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


/. k • k + 1 ' k 

*((£ V a kjj- 2 #(a i> 


2 0 


(B.7) 
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therefore: 


(B.8) 


a i - a k + l))= 2 / (a i 2) “^ (a k + l ) 


i ^f (a k + l a i ) 


g 0 


// k+ 1 \ 2 \ k+ 1 9 

iLh'jp'Zf'*' 


then, by mathematical induction: 


^(.f i (#)= 2 I / (a i 2) 


(B.9) 


(B. 1.0) 


QED. 
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appendix c 


CONVENTIONAL AUTOMATIC LANDING SYSTEM DESIGN 


C.l Introduction 


The most common type of automatic landing system uses information derived 
rom the IDS localizer and glide-path receivers to provide ‘the position error data 
required for acquisition and flight along the reference line defined by the intersection 
ox the ILS localizer and glide-slope center planes. Conventional systems provide 
essentially linear position control relative to this reference line. This appendix 
describes lateral and vertical control system designs which are based on the 
conventional approach. <. 

C ' 2 Conventional Lateral Position Control System 

In the lateral channel the coupler output provides the reference input to the 
roll autopilot. Roll angle serves as the lateral control variable, since the heading 
rate tjj is roughly proportional to roll angle <j> 


l = g tan t 


h 


(C.2-1) 


and the lateral velocity y is proportional to heading angle relative to the path (in 
the absence of a cross wind). 


y = VpSin^ 


(C.2-2) 


where g is the gravitational constant and v is the path velocity. A schematic diagram 
of a conventional LATERAL control system is shown in Fig. C.2. 1 where: 


K 


K. 


K. 

ty 


V 


loc 


is an adjustable position feedback gain. 

is an adjustable rate feedback gain. 

is the integral compensator gain. 

is the ILS Receiver Time Constant. 

is the velocity lead network time constant. 

is the estimated distance to the localized antenna. 
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Fig. C.2. 1 Conventional lateral position control system 







Position information relative to the localizer beam center is obtained by 
multiplying the angular deviation (radians) provided by the localizer receiver by 
the distance dj to the localizer antenna. This position information is multiplied 
by the gain K^to close the position control loop. 

The integral compensator operates on y to provide a roll command which 
maintains the correct crab angle when the aircraft is operated in a varying cross 
wind as well as correcting for any individual errors in trim. 

The lead network generates a signal proportional to the rate of change of 
lateral position; this signal provides dynamic response compensation as indicated 
above. This. compensation may be augmented by heading feedback due to the relation 
in Eq (C.2-2). 

The open-loop transfer characteristics of the control system are shown in 
Fig. C.2. 2, which uses the linearized vehicle transfer functions in Appendix B, ref 26, 
the roll angle control system in Appendix D and the parameter values in Table C.3-1. 
The closed- loop transfer function is shown in Fig. C.2. 3. A linearized response 
appears in Fig. C.2. 4. 

C. 3 Conventional Vertical Position Control System 

Vertical position control is achieved by applying the output of the vertical 
coupler as a reference input to the pitch autopilot as shown in Fig. C.3. 1 

where 

K z is an adjustable position feedback gain. 

is an adjustable velocity feedback gain. 

K. is an adjustable integral compensator gain, 
xz 

T^ is the ILS receiver time constant. 

T^ is the velocity lead network time constant. 

dg g is the estimated distance to the glide-slope antenna. 

a is the angular deviation from the glide-path center. 

z is the vertical coordinate of the glide-path center. 

g fa 

x is the distance to the glide-slope antenna. 

z is the vertical distance between the aircraft and the glide-path center. 

The vertical component of velocity z is approximately proportional to the pitch 0. 

z ~ -v 0 
P 
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Fig. C.2. 4 ' Unit-step response of the linearized conventional lateral 

position control system. 
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Fig. C.3-1 Conventional vertical position control system. 









The vertical distance between the aircraft and the glide-path center is obtained by 
multiplying the angular deviation (radians) from the ILS glide-slope receiver by the 
estimated distance to the glide-slope transmitter. This is then multiplied by a 
gain K z to close the position loop. 

The integral compensator provides the constant component of Q which is 
required for flight down an inclined path with zero position error as well as cor- 
recting for trim. The velocity lead network provides a signal proportional to the 
rate of change of vertical position for dynamic compensation. 

The open- loop transfer function of the control system is shown in Fig. C. 3. 2 
and is generated using the vehicle transfer function in Section C. 5, the pitch angle 
control system in Section 5. 11 and the parameter values in Table C.3-1, of ref 26. 
The closed-loop transfer function is shown in Fig. C. 3. 3, and the unit-step 
response of the linear model appears in Fig. C.3.4. 


Table C.3-1 Parameter Values for Conventional ILS Coupler 


GAINS 

K y 

lateral position gain 

0. 0205 deg/ft 

K- 

y 

lateral velocity gain 

0 , 410 deg/ft/ sec 

K. 

iy 

y intrgral compensator 
gain 

• 0. 005 sec 

K 

z 

vertical position gain 

0.041 deg/ft 

*z 

vertical velocity gain 

0. 205 deg/ft/sec 

K. 

1Z 

z integral compensator 
gain 

0. 005 sec 

CONSTANTS 

T 

r 

ILS receiver time 
constant 

0. 40 secs 

T v 

velocity filter time 
constant 

1. 00 secs 
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Fig. C. 3. 2 Open-loop transfer function of a conventional vertical 
position control system. 
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Fig. C.3. 3 Closed-loop transfer function of a conventional vertical 
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position control system. 
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Fig. C.3. 4 Unit-step response of the linearized conventional vertical 
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APPENDIX D 


ATTITUDE ’CONTROL SYSTEMS 


D. 1 Introduction 

As indicated in chapter 3 control of roll angle cp and pitch angle 0 plays an 
essential role in the construction of the lateral and vertical position control systems. 
A complete discussion of roll and pitch control is found in chapter 5 of reference 26, 
The roll and pitch control system designs are summarized here for completeness. 

D. 2 Roll Angle Control System Design 

A linearized model of the roll control system is shown in Fig. D. 2. 1. The 

open and closed loop magnitude and phase frequency domain characteristics are 

shown in Figs. D. 2. 2 and D. 2. 3 and the linearized system step response is shown 

in Fig. D. 2. 4. The complete roll control schematic in Fig. D. 2. 5 includes a lateral 

Stability Augmentation System (SAS) which generates the rudder deflections required 

for turn coordination. (0 = 0) Since the control system, with integral compensation 

(K. # 0), is conditionally stable and the roll control system is subject to saturation 
1CD 

the integral compensator output must be limited to ensure stable operation. The 
parameter values associated with Fig. D. 2. 4 are given in Table D.2-1 . 

D. 3 Pitch Angle Control System Design 

The pitch angle control system appears in linearized form in Fig. D. 3. 1. 
Phase and magnitude frequency domain characteristics of the linearized model are 
shown in Figs. D. 3. 2 and D. 3. 3. The step response of the linearized system is 
shown in Fig. D. 3. 4. The complete -pitch control system is shown in Fig D. 3. 5. 
The parameter values associated with Fig. D. 3. 5 are shown in table D. 2-1. 
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Fig. D. 2. 1 Linear model of the roll angle control system. 
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Fig. D. 3. 2 Open-loop transfer function of the pitch control system. 
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Fig. D.3. 5 Pitch attitude control system. 
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APPENDIX E 

AIRSPEED CONTROL SYSTEM 


E. 1 Introduction 

The basic nonlinear vehicle equations were derived in chapter 8. Linear- 
ization about the nominal operating condition provides for application of the tech- 
niques of linear control system analysis and synthesis. 

E. 2 Derivation of the Airspeed-Thrust Transfer Eunction 

The longitudinal variables u, -a, q are related by a coupled set of linear 
differential equations in the perturbation quantities as follows: 


C-u = C*u-i-C*ar + C -6 + C*. 6„ 

u uu - uor ■ u0 u5 rp T 


(E. 2-1) 


C • a = C . u + C* a 4- C« q + CmS 
a cru a a <*q a-0 


-i-C*,. 5 + C 5rp 

e T 


(E. 2-2) 


C • q = C • u 4- C • 0 " + C* q + C*. 6 

q M qu q« qq qs e e 


+ C * *a 4- C* 

qa- q6 T T 


{E. 2-3) 


where, for example, is the coefficient of the u. perturbation in the u equation. 
Airspeed control is achieved by throttle manipulation; the resulting thrust variation 
couples into the pitch and lift equations. To obtain an approximate transfer function 
relating u and 6^., the longitudinal equations are simplified for the moment by 
assuming that the pitch control system maintains a zero pitch rate through elevator 
manipulation. 

The simplified equations are: 


Cm 

u 


C- u 
uu 


4- C • a 
u a 


+ C 


U6 r 


(E. 2 4). 


C • a ~ C. u 4- C • » 4- C *, S,-, 4- C.j. 5 

a cm aa T a5 e 

i e 


(E. 2 5) 
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(E. 2-6) 


0 = C. u + C. a + C..a + C.. 5„. + C. 6 

qu qa qa q6 T T r q« e e 

These equations may be rewritten in Laplace notation in the form 


G u + G a 
uu ua; 


= G u 5 t - 6 T 


G u + G # + G <, 6 = G . 6™ 

au act a8 e e 0-6 T 


G u + G a + G . S„ = G_, 5 


4 « I U 1 r V 

qu q« q6 e 


q« T T 


(E. 2-7) 
(E. 2-8) 
(E. 2-9) 


where the G uu , G uqi .... are polynomials in s. 
Defining the characteristic polynomial 
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permits a simultaneous solution of (E. 2-7, 8, 9) for the transfer function 


!£U ’ 


G q6„ G q«]- G u»[ G «6 T G q6 ' G * 6 G qi J 
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E. 3 Numerical Value of Airspeed -Thrust Transfer Function 

The numerical value for the transfer function developed in the previous 
section incorporating coefficient values appropriate to the Boeing B2 707- 200 is 


ru 1 _ 0.0000824 (s + 0.668) 

l^T Jq =0 (s + 0.653) (s +0.053) 


E. 4 Airspeed Control System Design 

The design of the airspeed control system is shown in block diagram form 
in Fig. E. 4. 1. Airspeed and a lagged acceleration are the feedback signals supplied. 
The first blocks in the forward loop comprise a proportional plus integral network 
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to provide compensation. A first order filter is provided to limit high frequency 
throttle activity and a first order equation models the turbine engine thrust lag. 

The engine time constant was assumed to be 4 seconds. The time constants in the 
filter and the acceleration feedback loop were each set at 1/10 second. 

The root locus for the inner (acceleration) loop is shown in Fig. E. 4. 2. The 
closed loop poles were chosen on the real axis, as shown, and the corresponding 
value for the gain was 


K* = 10,200 

u 

With this value for K« and with K. « 0 the root locus for the outer loop is 

u xu ^ 

shown in Fig. E. 4. 3. The dominant closed loop poles were selected to give a 
damping ratio of about £ =0.7, The corresponding value for was 

K = 5, 650 
u 

The control system was incorporated with the above gains into the digital 
simulation, which includes the full nonlinear aircraft equations of motion, and 
subsequently no changes were found to be necessary. The simulation also deter- 
mined that a satisfactory compromise value for the P +1 compensation gain was 

K. = 0. 02 

xu 


E. 5 Airspeed Control System Transfer Functions 
Define 


r 1 i 

r 1 i 

no 

Lt.s + Ij 

_T S + 1 J 

-O 


H = 


r s +K. r K*s -i x 

K [-5^“] + fern] } 


u 


n r- S 4- IC. -i 

G [-s-“K 


1 +GH 


For the design gains, a computer program has factorized the above transfer 
functions and developed the corresponding Bode plots. 





Fig. E 4. 3 Root locus of the airspeed control system with K^ = 10, 200 
and K^ u = 0. (Root loci not to scale) 


E-5 





0.000206 (s + 0.668) 

" (s + 10) (s +0.25) (s + 0.653) (s + 0.053) 


107,650 (s + 0.02) (s + 0.5) 
s(s + 10) 


' -u - 1.162 (s + 0.02) (s + 0.67) (s + 10) 

u d <s + 0.0194) (s + 0.253 + 0. 2531) (s + 0.65) (s + 8.40) (s +11.38) 


The vehicle transfer function, G, is plotted in Fig. E. 5. 1, the feedback 
function, H, in Fig. E. 5,.2, and the closed loop transfer function, u/u^, in Fig. 

E. 5. 3. The unit step response of the transfer function u/u^ is shown in Fig. E.5. 4. 
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Fig. E.5. 3 Bode plot of transfer function, u/u 
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